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I.  Program  Overview 


I 
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This  final  report  contains  the  first  year's  results  of  a project 
originally  proposed  as  a three-year  research  program.  The  main  objec- 
tive was  the  theoretical  and  experimental  investigation  of  the  propagation 
and  scattering  of  ultrasonic  waves  in  fiber-reinforced  composite  materials, 
i both  with  and  without  flaws.  The  program  was  supported  by  AFSOR  from 

April  1,  1976  to  March  31,  1977. 

It  was  recognized  at  the  outset  that  composite  materials  present 
unique  ultrasonic  inspection  problems.  First,  there  is  their  inherent 
inhomogeneity.  For  unidirectional  fibers  the  propagation  of  elastic  waves 
^ will  be  similar  to  that  in  anisotropic  materials  such  as  single  crystals. 

Second,  there  is  a characteristic  dimension,  the  filament  diameter  or 
its  ply  spacing,  associated  with  the  microstructural  aspects  of  the  com- 
posite. For  a boron-epoxy  composite  these  are  0.08  mm  and  0.20  mm,  respec- 
^ tively.  Since  ultrasonic  poises  used  in  most  NOT  situations  contain  fre- 

quency components  in  the  megahertz  range,  interrogation  pulses  will  con- 
tain wavelengths  on  the  order  of  the  material  dimension,  and  the  propaga- 
tion of  ultrasonic  pulses  through  the  composite  will  be  greatly  affected 
by  the  multiple  scattering  of  waves.  The  observed  dispersion  and  attenua- 
tion of  ultrasonic  waves  in  composites  is  a result  of  the  multiple  scat- 
tering. Finally,  in  many  cases,  particularly  when  the  composite  consists 
of  an  epoxy  matrix,  the  ultrasonic  signals  are  highly  damped. 

The  research  program  completed  has  focused  on  each  of  the  above  areas, 
and  we  believe  that  significant  progress  has  been  made.  The  research  pro- 
gram has  benefitted  from  close  interaction  between  the  experimental  and 
the  theoretical  work,  and  from  support  received  from  other  sources  prior 
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to  the  initiation  of  this  program.  Experimental  techniques  utilizing 
continuous  waves  and  broadband  pulse  measurements  were  developed.  Analyses 
of  the  multiple  scattering  of  elastic  waves  by  fibers  imebdded  in  other-  ' 

wise  homogeneous  matrix  materials  were  made.  Details  are  described  in  ' 

the  next  section. 

During  the  year,  two  faculty  members  (W.  Sachse  and  Y.H.  Pao),  two  | 

post-doctoral  research  associates  (V.V.  Varadan  and  C.S.  Ting)  and  sev-  I 

i 

eral  graduate  students  worked  either  part-  or  full-time  on  this  project. 

II.  ACCOMPLISHMENTS  j 

(A)  The  frequency-dependent  propagational  characteristics  of  ultra-  | 

sonic  waves  in  various  directions  of  boron-spoxy  specimens  have  been  in- 
vestigated. Two  ultrasonic  techniques  were  developed  for  the  accurate  i 

t 

determination  of  ultrasonic  phase  and  group  velocities  in  highly  absorb-  j 

tive  composite  materials.  Either  technique  can  also  be  used  for  the  deter-  ! 

mination  of  the  dispersion  relation  for  material.  Both  techniques  have  been  | 

demonstrated  with  longitudinal  and  shear  waves  of  frequencies  ranging  from  | 

I 

j 

0.5  to  about  10  MHz  propagating  in  various  directions  to  the  fiber  axis.  | 

r 

i 

The  continuous-waves  technique  is  an  improvement  on  the  existing  j' 

method  of  phase  comparison.  By  transmitting  with  variable  frequency,  one  ^ 

can  determine  the  phase  velocity  of  waves  in  a composite  if  the 
number  of  complete  sinusoidal  waves  is  known.  In  this  method,  a constraint 
equation  is  found  which  permits  the  determination  of  the  absolute  number 
of  wave  lengths  of  a wave  in  a specimen. 

This  work  was  presented  at  the  89th  Meeting  of  the  Acoustical  Society 
of  America  at  the  Pennsylvania  State  University  in  June  1977  and  a paper 
(copy  appended)  has  been  submitted  for  publication  to  the  Journal  of  the 
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Acoustical  Society  of  America. 

(B)  The  second  technique,  ultrasonic  phase  spectroscopy,  measures 
the  phase  spectrum  of  a broadband  ultrasonic  pulse.  We  have  shown  that 
the  phase  spectrum  of  a propagating  pulse  is  related  to  the  dispersion 
relation  of  the  medium  from  which  the  phase  and  group  velocities  of  the 
medium  can  be  calculated  as  a function  of  fequency.  A digital  phase 
spectrum  analyzer  has  been  set  in  operation  to  determine  the  propagational 
characteristics  of  longitudinal  and  shear  broadband  ultrasonic  pulses  pro- 
pagating in  boron-epoxy  specimens  and  in  other  dispersive  materials. 

This  portion  of  the  program  was  supported  in  part  by  the  National 
Science  Foundation  through  grants  to  the  College  of  Engineering  and  the 
Materials  Science  Center  at  Cornell  University.  The  phase  spectroscopy 
technique  was  described  in  a presentation  at  the  ARPA/AFML  Review  of  Quanti- 
tative NDE  in  June  1977.  A paper  (copy  appended)  has  been  accepted  for 
publication  in  the  Journal  of  Applied  Physics. 

(C)  In  the  theoretical  study  of  the  multiple  scattering  of  waves  by 
fibers  in  a composite,  the  newly-developed  method  of  the  transition  matrix 
was  extended  to  the  case  of  multiple  scattering  of  Sh  waves  (scalar  waves) 
normally  incident  on  cylinders  of  arbitrary  cross  section.  Using  the 
transition  matrix  of  a single  obstacle,  statistical  averaging,  and  Lax's 
"quasi -crysstal 1 ine"  approximation,  we  developed  equations  for  the  average 
amplitudes  of  scattered  and  excitation  fields,  from  which  the  dispersion 
relation  and  hence  the  phase  and  group  velocities  for  the  material  are 
obtained.  In  the  Rayleigh  (low  frequency)  limit,  the  effective  elastic 
properties  of  the  composite  agree  with  the  existing  results  in  the  literature. 

These  results  were  reported  in  the  Compoiste  Materials  Review,  held  at 
Wright-Patterson  AFB  in  October  1976;  and  they  have  been  submitted  for 


III.  CONCLUSIONS 

The  propagation  and  scattering  of  ultrasonic  waves  and  pulses  in 
composite  materials  have  been  studied.  Two  experimental  techniques  have 
been  developed  for  accurately  determining  the  dispersive  nature  of  the 
propagation  of  ultrasonic  waves  in  such  materials.  A theory  for  the 
multiple  scattering  of  elastic  waves  by  thin  fibers  of  arbitrary  shape  was 
formulated.  Based  on  the  theory,  dispersion  relations  are  found  for 
the  ultrasonic  longitudinal  and  shear  waves  propagating  normal  to  the 
fiber  directions. 
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MEASUREMENT  OF  ULTRASONIC  DISPERSION  BY  PHASE  COMPARISOiJ 


OF  CONTINUOUS  HARMONIC  WAVES 

C.  S.  Ting  and  Wolfgang  Sachse 

Department  of  Theoretical  and  Applied  Mechanics 
Cornell  University,  Ithaca,  New  York  - 14853 

Abstract 

The  method  of  phase  comparison  of  continuous  waves  is  applied  to 
determine  the  dispersion  relation,  phase  and  group  velocities  as  a 
function  of  frequency  of  dispersive  materials.  A combination  of  the 
variable  frequency  method  and  the  variable  path-length  method  is  found 
necessary  to  eliminate  any  uncertainty  in  the  dispersion  relation 
determination.  Experiments  are  performed  on  specimens  of  various 
thicknesses.  A constraint  equation  can  be  derived  since  the  dispersion 
relation  is  unique  and  independent  of  the  specimen  thickness.  This 
equation  provides  a procedure  for  determining  the  absolute  number  of 
wavelengths  in  the  specimen.  A transducer  and  electronics  compensation 
procedure  is  also  used.  Measurements  in  uni-directional,  fiber-reinforced 
Boron-Epoxy  specimens  show  good  agreement  with  results  reported  previously. 
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We  describe  in  this  paper  a continuous-wave,  through-transmission, 
phase  comparison  method  for  measuring  the  phase  and  group  velocities  of 
ultrasonic  waves  in  dispersive  materials.  The  method  is  an  adaptation 
of  the  phase-comparison  methods  described  by  Lynnworth,  et  aj_  [1]  which 
have  been  reviewed  by  Papadakis  [2].  The  procedure  we  shall  describe 
provides  an  unambiguous  means  for  determining  the  exact  wavelength  of 
an  ultrasonic  wave  in  a solid  specimen.  This  is  then  used  to  find  the 
dispersion  relation,  phase  and  group  velocities  as  a function  of  frequency 
for  the  material. 

I.  Introduction 

The  speed  of  propagation  of  a monochromatic  wave  can  be  defined  as  the 
speed  with  which  the  phase  of  the  wave  is  propagated  [3].  A plane  wave, 
propagating  undamped  in  the  x-direction,  can  be  described  at  (x,t)  as 

u(x,t)  = (1) 

where  the  quantity  (oittkx)  is  known  as  the  phase  function,  a)(=2:rrf)  is 
the  circular  frequency  and  k (=2Tr/X)  is  the  wavenumber.  Surfaces  of 
constant  phase  propagate  at  speed  v where 

V = dx/dt  = u)/k  = fX  (2) 

V is  the  phase  velocity. 

It  can  be  shown  [4]  that  a pulse,  being  a superposition  of  waves, 
has  its  mean  amplitude  propagate  with  a group  velocity  given  by 


U = dw/dk  = df/d(l/X) 


(3) 


11  wtm 
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When  a wave  is  subjected  to  frequency-dependent  damping,  complications 
arise  in  defining  its  phase  and  group  velocity.  These  are  discussed  by 
Sachse  and  Pao  [5]. 


From  Equations  (2)  and  (3)  one  obtains 


U = 

^ dk 


V ^ 

dk 


(4) 


an  equation  relating  the  phase  and  group  velocities.  Thus,  only  when 
V f v(k)  will  the  group  velocity  equal  the  phase  velocity.  When  the 
propagating  medium  is  dispersive,  the  phase  and  group  velocities  are 
unequal  and  both  depend  on  the  frequency  (wavelength)  of  the  wave  Further- 
more, pulses  change  their  shape  while  propagating  through  such  m'  jia. 


The  measurement  of  phase  and  group  velocities  with  continuous  waves 
and  pulse  techniques  in  dispersive  materials  has  been  reviewed  by  Young 
[6],  Papadakis  [2J  and  recently  by  Sachse  and  Pao  [5].  This  paper  is 
restricted  to  a discussion  of  measurements  made  with  continuous  harmonic 
waves  in  dispersive  solids. 


The  phase  velocity-frequency  behavior  and  hence  the  dispersion  relation 
can  be  determined  from  continuous  wave  measurements  only  if  the  wavelength 
of  the  wave  is  known  unambiguously.  This  is  usually  impossible  in  solids 
unless  they  are  transparent  so  that  the  wavelength  of  a wave  can  be  directly 
measured  as  a function  of  frequency  with  optical  techniques  [7]. 


In  an  attempt  to  measure  the  wavelength  of  the  wave,  Lynnworth,  et  al 
[IJ  used  specimens  consisting  of  two  wedges  such  that  the  propagation  path 
length  between  source  and  receiving  transducer  could  be  varied.  The  lengtn 
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changes  corresponding  to  a 2it  phase  change  in  the  received  signal  relative 
to  the  input  signal,  should  correspond  to  the  wavelength  in  the  specimen. 

The  difficulty  with  the  technique  appears  to  be  associated  with  anomalous 
effects  arising  from  the  interface  and  microstructure  of  the  materials  used. 

The  group  velocity  can  be  determined  once  the  dispersion  relation 
k = k(aj)  (or  d)  = u)(k))  or  the  frequency  dependence  of  the  phase  velocity  is 
known  for  the  material.  If,  however,  only  the  group  velocity  is  known  as  a 
function  of  frequency,  it  is  impossible  to  determine  v uniquely  by  integra- 
tion of  Eqn.  (4).  Missing  is  the  precise  value  of  k at  the  lower  limit  of 
integration. 

Because  a wavelength  determination  solids  is  difficult,  in  a continuous 
wave  experiment  the  frequency  of  excitation  is  varied  and  the  resonance 
frequencies  of  the  specimen  determined.  Even  though  the  technique  has  been 
in  existence  for  a long  time,  few  measurements  in  dispersive  media  can  be 
found  in  the  literature.  One  example  previously  cited,  is  the  work  of 
Lynnworth,  £t  ^ [1]  who  measured  the  group  velocity  this  way  in  3-dimensional 
carbon-phenolic  specimens. 

The  measurement  technique  is  sometimes  called  the  7r-phase  comparison 
technique.  In  it,  a continuous  wave  of  freqnency  f is  transmitted  through 
a specimen  of  length  L.  Data  are  obtained  by  determining  the  lowest 
frequency  at  which  the  signal  received  at  L is  out-of-phase  with  the 
input  signal.  This  corresponds  to  a half-wavelength  in  the  specimen,  hence 
A/2  = L.  Increasing  the  frequency  of  excitation,  the  received  signal  is 
successively  brought  in-phase  and  Tr-radians  out-of-phase  with  the  input 
signal.  At  each  condition  the  number  of  cycles  of  the  wave  is  increased 
by  one-half  and  so  the  number  of  cycles  in  the  specimen,  N,  is  known  at 


every  frequency.  The  process  is  continued  until  no  signal  is  detected 
at  the  transducer,  being  the  result  of  specimen  attenuation  effects  or 
transducer  bandwidth  limitations. 

The  dispersion  relation  and  the  phase  velocity  are  determined  at 
each  frequency  f,  for  which  N is  known.  Since 

X = L/N  (5) 

the  dispersion  relation  k or  oj  is  specified.  The  phase  velocity  is 

V = fX  = (6) 

In  addition,  the  group  velocity,  U , is  known,  since 

U = df/d(l)  = L df/dN  (7) 

The  crucial  point  in  the  wavelength  determination  is  that  N is  known 
unambiguously  over  the  entire  frequency  range.  If  this  is  not  the  case, 
the  group  velocity  U can  still  be  determined  as  follows.  If  at  frequency 
f^  an  in-phase  or  an  out-of-phase  condition  exists  such  that  N = N-j  and 

if  the  frequency  is  then  adjusted  to  f2  where  ^2  ” ^1  ^ 
which  another  in-phase  or  out-of-phase  condition  exists  such  that 
1^2  = + AN,  then  we  find 


Ak  ■ k2-k^  ~ ^ 27rAN 

In  the  limit,  as  k 0,  this  equals  the  group  velocity  U.  Thus, 

U = L df/dN  (9) 


It  was  pointed  out  in  Ref  (6)  that  it  would  be  possible  to  determine 


the  phase  velocity  from  Eqn.  (4)  by  integration,  provided  that  the 
2irN, 

value  for  k^(=  — j— which  appears  in  the  integration  constant  is 
known.  This  is  in  agreement  with  the  discussion  earlier. 

A unique  determination  of  N,  the  number  of  wavelengths  in  a speci- 
ment,  is  experimentally  quite  difficult  if  not  impossible  in  non-transparent 
I solids.  As  Papadakis  has  pointed  out  [2],  at  low  frequencies  spurious 

interferences,  sidewall  reflections  and  complex  modes  in  a specimen  make 
it  difficult  to  identify  the  first  cut-of-phase  condition  as  well  as  sub- 
sequent in-  and  out-of-phase  conditions.  With  composite  material  specimens, 
we  have  observed  the  multimode  phenomena  to  occur  below  about  0.3  MHz  which 
^ appears  to  be  inherent  to  such  materials.  The  presence  of  the  phenomena 

is  indicated  by  several  successive  in-phase  (or  out-of-phase)  conditions 
oeing  observed  at  irregular  frequency  intervals. 

If  the  material  is  non-dispersive  in  some  frequency  range,  then  v = U 
^ ana  a group  velocity  measurement  can  be  used  to  determine  the  wavelength 

and  hence  H at  a frequency  in  that  range.  If,  however,  the  medium  posses- 
ses stop-bands  [8],  the  wavelength  becomes  indeterminate.  In  the  next  sec- 
tion we  describe  a method  which  combines  the  variable  frequency  method 
and  the  variable  path-length  method  to  determine  the  exact  number  of 
cycles,  N,  in  a specimen  as  a function  of  frequency. 

II . Basis  of  the  Method 

The  measurements  on  a specimen  are  begun  at  any  frequency  f at  which 

the  input  and  output  signals  are  in-phase.  This  condition  is  assigned  any 

arbitrary  integer  value.  The  excitation  frequency  is  then  increased  and 

★ 

the  number  of  cylces  of  phase,  N , relative  to  the  initial  point,  recorded 


as  a function  of  frequency.  It  is  clear,  that  in  order  to  obtain  the 
absolute  value  of  cycles  of  phase,  an  integer  value  of  variation,  6N,  must 


J 
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be  added  to  each  point. 

The  number  of  relative  cycles  of  phase,  however,  also  contains  the 

effects  of  frequency-dependent  phase  shifts  caused  by  the  transducers 

and  amplifiers.  This  is  denoted  by  n(f)  and  it  must  be  subtracted  from 
★ 

the  N + 6N  cycles  of  relative  phase  to  obtain  the  absolute  number  of 
cycles  of  phase  in  the  specimen.  Hence 

N = N*  + 6N  - n(f)  (10) 


A.  n(f)  Determination 

If  the  delay  time  t^  resulting  from  the  transducers  and  amplifiers 
is  a constant  as  a function  of  frequency,  then  the  cycles  of  phase 
shift  resulting  from  these  is 


n 


(11) 


More  generally,  however,  the  cycles  of  phase  shift  caused  by  these  elements 
is  also  frequency-dependent,  hence  n(f)  must  be  determined  experimentally 
Experience  has  shown  that  determination  of  n{f)  is  to  be  preferred  over 
the  use  of  a second  set  of  transducers  which  are  in' intimate  contact  and 
coupled  to  an  amplifier.  It  is  generally  very  difficult  to  obtain  two 
sets  of  ultrasonic  transducers  whose  frequency-dependent  phase  shifts  are 
identical . 

The  determination  of  n(f)  can  be  accomplished  in  two  ways.  In  one, 

the  transducers  are  placed  in  intimate  contact  and  the  function  n(f)  is 

★ 

determined  directly  from  the  cycles  of  phase  shift,  n , measured  between 
input  and  output  transducers  as  a function  of  frequency.  Alternately,  a 
non-dispersive  highly-attenuating  specimen  is  used  in  the  phase-compari- 
son system.  This  facilitates  measurements  to  be  reliably  made  to  the 


lowest  out-of-phase  frequency.  From  Eqn.  (5) 


N = L/X  = ^ 

If  the  phase  velocity  of  the  non-dispersi ve  material  is  v^,  then  the 
aispersion  relation  is 


k = = f=-f  (12) 

0 0 

Combining  these  equations,  one  finds 

N = (13) 

0 

Therefore,  the  non-specimen  phase  shift  is  determined  from 

n(f)  - n*(f)'-^f  (14) 

0 

and  this  must  be  subtracted  from  the  measured  cycles  of  relative  phase 
Shi  ft. 

B.  6N  Determination 

The  value  of  5N  can  be  determined  by  performing  phase-comparison 

measurements  on  specimens  of  various  lengths.  Since  the  wavenumber  is 

directly  related  to  the  number  of  cycles  of  phase  in  a specimen  according 

to  Eqn.  (5),  then,  if  a set  of  m specimens  of  length  is  tested  and 
★ 

for  each,  N^(f)  is  measured,  the  following  set  of  m equations  result 

ki(f)  = ^N*(f)+5N.-n(f)]  (15) 

The  dispersion  relation,  w = u)(k),  is  unique  for  the  material  and  independent 
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of  specimen  length.  It  follows,  therefore,  that  at  any  frequency  f 

k^Cf)  = k2(f)  = k3(f)  * ...  = k^(f)  (16) 
or 

^N*(f)+6N^-n(f)]  = ^iN*(f)+6N2-n(f)]  = ...  = ^N*(f)+6N^-n(f )]  (17) 

1 2 ^m 

Equations  (17)  are  quivalent  to  a set  of  (m-1)  linearly  independent 
equations  for  the  m unknowns,  6N, , 6No,  6N-,  ...  5N  . An  additional 
condition  is  that  the  6N^  are  to  be  integer  values.  Hence,  the 
can  all  be  uniquely  detenrnned. 

Eqns.  (17)  can  be  rewritten  in  the  form 

6N.  = 6N,  + B.  i = 2,  3,  4,  ...,  m (18) 

1 L ^ I I 

where 

= {^N*(f).L^N*(f)-(L.-L^)n(f)} 

For  6N^  to  be  independent  of  frequency,  requires  that  be  also.  The 
requirement  that  the  values  be  exact  integers  is  generally  not  attain- 

able frdm  experimental  measurements.  The  requirement  is  therefore  relaxed 
such  that  the  6N^  be  close  to  integer  values.  Thus,  in  practice,  one 
finds  the  set  of  6N^  which  satisfy  Eqn.  (18)  and  whose  summed  deviation 
of  each  5N^.  to  its  nearest  integer  is  a minimum. 

The  solutions  obtained, depend  clearly  on  the  various  specimen  lengths 
chosen.  In  order  to  obtain  a solution  without  ambiguity,  one  must  chose 
specimen  lengths  which  preclude  multiple  solutions. 

An  extreme  case  is  that  in  which  all  the  specimens  are  of  equal  length. 
In  that  case,  B^.  = 0 and  5N.  cannot  be  determined  uniquely.  The  optimal 
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choice  of  specimen  lengths  to  be  used  in  a series  of  experiments  is  not 
immediately  apparent.  If  the  thickness  ratios  between  specimens  is 

....  the  common  denominator  is  as  large  as  possible.  If  five 
specimens  are  considered,  with  * 1 , and  if  at  a value  of  6N^  the 
other  are  all  integer  values,  only  when  SN-j  is  equal  to  60  will  all 
the  other  equal  integer  values  again.  It  is  clear  that  a judicious 

choice  of  specimen  thickness  ratios  should  be  used  if  the  dispersion  rela- 
tion for  the  material  is  to  be  determined  unambiguously. 

Although  the  number  of  specimens  is  at  minimum  two,  a larger  number 
of  specimens  greatly  simplifies  the  determination. 

Ill . Experiments 

The  electronic  system  for  making  ultrasonic  phase-comparison  measure- 
ments is  shown  in  Figure  1.  It  is  similar  to  that  described  by  Papadakis 
[2]  and  consists  essentially  of  a continuous  radio-frequency  oscillator 
(HP  606A),  a frequency  counter  (HP5326B),  transducers  (Panametrics  VIP- 
Series),  specimen  holder,  amplifiers  and  display  oscilloscope  (Tektronix 
7A18,  7704A).  Both  input  and  received  signals  are  displayed  on  the  oscil- 
loscope. The  signals  were  amplified  to  the  same  level  and  added.  A null 
signal  was  used  to  indicate  an  out-of-phase  condition.  The  in-phase  condi- 
tion was  determined  similarly  with  one  of  the  signals  inverted.  A vector 
voltmeter  can  also  be  connected  to  input  and  receiving  transducers  and 
the  phase  measured  directly. 

The  determination  of  the  tranducer  and  electronic  phase  delay,  n(f), 
was  determined  by  the  two  methods  described  in  Section  II.  The  results  are 
shown  in  Figure  2.  The  data  points  indicated  by  Method  (I)  are  the  results 
obtained  when  the  two  transducers  were  in  intimate  contact  and  the  phase 
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of  the  received  signal  was  measured  relative  to  the  excitation  signal  as 
a function  of  frequency.  The  line  indicated  on  Figure  2 passes  through  the 
points  determined  from  the  data  of  a 1.27  cm  thick  specimen  of  plexiglass 
placed  between  the  transducers  with  the  data  being  reouced  according 
to  the  procedure  described  in  Section  II.  The  measurements  were  repeated 
with  a 1.27  cm  thick  specimen  of  brass.  The  results  were  identical. 

To  demonstrate  the  applicability  of  the  data  analysis  procedure, 
phase-comparison  measurements  were  made  on  specimens  cut  form  a 96-ply 
"Scotchply"  uni-directional,  boron-epoxy  composite  panel.  The  specimen 
material  is  the  same  as  that  used  in  some  earlier  elastic  moduli  experi- 
ments and  the  characteristics  of  the  material  have  been  reported  previ- 
ously [9].  Shear  waves  were  used  in  the  experiments  to  be  described. 

The  wave  propagation  direction  was  coincident  with  the  fiber  direction 
and  the  particle  polarization  direction  was  perpendicular  to  the  ply  layers 
of  the  composite. 

Five  specimens  of  the  same  orientation,  but  possessing  different 

thicknesses  were  used.  The  original  data  obtained  from  phase-comparison 

* 

measurements  is  shown  in  Figure  3.  The  number  of  cycles  of  phase,  N , 
indicated  on  the  ordinate  is  measured  relative  to  a particular  in-phase, 
out-of-phase  condition.  The  values  of  in  Eqn.  (18)  were  evaluated 
at  five  frequencies:  2,  3,  4,  5,  and  6 MHz.  The  values  of  N.  at  each 
frequency  were  obtained  by  a linear  interpolation  of  the  two  adjacent 
data  points.  With  the  values  of  n(f),  it  was  found  that  the  variation 
of  at  these  five  frequencies  is  less  than  0.07  cycle  of  phase  delay. 
Using  the  averaged  values  of  B^  and  the  length  ratios,  Eqn.  (18) 

is  used  to  evaluate  the  various  (i  > 1)  in  terms  of  the  6N^ . One 
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finds 

(SNg  = 0.911  + 0.563  (±0.035) 

5N3  = 0.310  6H^  + 0.396  (±0.021) 

(19) 

6N4  = 0.501  + 0.955  (±0.023) 

6N5  = 1.466  5N^  + 0.820  (±0.034) 

This  set  of  equations  is  plotted  in  Figure  4. 

As  stated  in  Section  II,  in  order  to  determine  the  appropriate  value 
for  the  5N.|,  the  sum  of  the  deviation  of  each  6N^. (i  > 1)  from  its  nearest 
integer  value  must  be  a minimum.  This  can  be  done  numerically  using  Eqns. 
(19)  or  it  can  be  determined  directly  from  Figure  4 with  a calibrated 
scale.  The  deviations  of  each  (i  >1)  from  its  nearest  integer  value 
are  summed  for  every  integer  value  of  61!^  and  this  total  deviation  is 
plotted  versus  6N-j  as  shown  in  Figure  5.  It  is  clear  that  for  (SN.|  = 2 
the  total  deviation  is  a minimum.  It  follows  then  that  5N2  ~ = 1, 

6N^  = 2 and  6Ng  = 4 for  the  example  considered. 

Once  the  have  been  determined,  the  dispersion  relation  for  the 

material  is  established  according  to  Eqn.  (15).  This  is  shown  in  Figure 
6 with  all  the  original  data  (Fig.  2)  now  fitting  onto  one  continuous  curve. 

The  phase  velocity,  v(=2TTf/k)  is  determined  directly  from  the  disper- 
sion curve.  The  result  for  the  0.543  cm  thick  specimen  is  shown  in  Figure 
7.  The  group  velocity  is  determined  by  differentiating  the  dispersion 
curve.  In  order  that  artificats  arising  from  numerical  differentiation 
of  the  discrete  points  of  the  dispersion  curve  be  minimized,  a polynomial 
is  fit  to  the  dispersion  data.  In  the  present  experiments,  a cubic  poly- 
nomial of  f as  a function  of  k was  used  with  the  resulting  polynomial 
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coefficients  given  in  Table  I.  The  group  velocity  is  then  determined  as 
a function  of  frequency  in  closed  form  from  the  dispersion  curve  polynomial. 
As  a check  of  this  procedure,  the  data  points  comprising  the  dispersion 
relation  were  analyzed  point-by-point  with  a four-point  Lagrangian 
interpolation  technique  DO]  to  determine  the  group  velocity.  The  results 
of  both  group  velocity  determinations  for  the  0.543  cm  thick  specimen  are 
also  shown  in  Figure  7. 

The  phase  and  group  velocity  determinations  made  on  all  five  speci- 
mens is  compared  in  Figure  8.  The  phase  velocity  results  are  nearly 
identical  while  the  group  velocity  results  show  a slightly  large  scatter. 

For  comparison,  the  results  of  Tauchert  and  Guzelsu  [11]  who  used  an  r.f. 
burst  measurement  technique  to  determine  the  time  delay  of  the  burst  and 
so  its  group  velocity  through  specimens  of  similar  material.  Although 
their  measurements  were  made  with  different  transucers  at  frequencies  of 
0.5,  2.25,  5.0  and  10  MHz,  the  line  shown  on  Figure  8 is  drawn  through  all 
their  data  points.  The  agreement  with  the  results  of  continuous  wave 
measurements  is  within  8%. 

IV.  Conclusions 

By  combining  the  variable  frequency  and  the  variable  path  length 
continuous  wave  techniques,  we  have  shown  that  the  dispersion  relation 
for  a material  can  be  determined  from  ir-phase  comparison  measurements. 

The  procedure  utilizes  phase  data  relative  to  an  arbitrary  starting 
frequency  which  has  been  modified  to  take  into  account  the  phase  delays 
associated  with  the  transducers  and  electronics.  Since  the  dispersion 
relation  for  a material  is  unique  and  the  number  of  cycles  must  be  an 
integer  number,  it  is  shown  that  the  number  of  cycles  at  the  test  start 
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can  be  established  and  so  the  dispersion  relation  for  the  material 
is  found  to  which  a polynomial  can  be  fitted.  Both  phase  and  group  veloci- 
ties can  then  be  determined  from  the  dispersion  relation  point-by-point 
numerically  or  in  closed  form  by  differentiation  of  the  polynomial. 

The  results  obtained  are  consistent  from  specimen  to  specimen, regard- 
. less  of  length.  Furthermore,  the  group  velocity  determination  is  in  agree- 

ment with  that  obtained  with  a r.f.  burst  technique. 


I 

I 


Acknowledgments 


This  work  was  supported  in  part  by  the  Air  Force  Office  of  Scientific 
Research  and  by  the  National  Science  Foundation  through  the  Materials 
Science  Center  at  Cornell  University. 


! 


li^ijipUHiyL^ 


References: 

1.  L.  Lynnworth,  E.  P.  Papadakis  and  W.  Rea,  "Ultrasonic  Measurement  of 
'Phase  and  Group  Velocity  Using  Continuous  Wave  Transmission  Techniques," 

AMMRC  Report  CTR  73-2j January  1973. 

2.  E.  P.  Papadakis,  Chapt.  5 in  Physical  Acoustics,  Vol . XII,  W.  P.  Mason 
and  R.  N.  Thurston,  Eds.,  Academic  Press, ^ew  York,  1976,  pp.  277-374. 

3.  R.  B.  Lindsay,  Mechanical  Radiation,  McGraw-Hill  Book  Co.,  New  York 
(1960),  ^ 27. 

4.  L.  Brillouin,  Wave  Propagation  and  Group  Velocity,  Academic  Press,  New 
York,  1960.  . 

5.  W.  Sachse  and#.  H.  Pao,  Report  #2773,  Materials  Science  Center,  Cornell 
University,  Ithaca,  NY  (July  1977).  To  appear  in  J.  Appl . Phys. 

6.  E.  H.  Young,  Jr.,^IRE  Trans.  UE-9,  13-21  (1962). 

7.  C.  H.  Palmer,  R.  0.  Glaus  and  S.  E.  Pick,  Applied  Optics,  1_6,  1849-1856 
(1977). 

8.  H.  J.  Sutherland  and  R.  Lingle,  J.  Composite  Matls.,  490-502  (1972). 

9.  W.  Sachse,  J.  Comp.  Matls.,  8,  378-390  (1974). 

10.  R.  W.  Hamming,  Numerical  Methods  for  Scientists  and  Engineers,  McGraw- 
Hill  Book  Co.,  New  York,  1962,  pp.  94-97. 

11.  T.  Tauchert  and  A.  N.  Guzelsu,  Trans.  ASME-J.  Appl.  Mech.,  98-102, 
(1972). 


# 


r 


- A - 15  - 


I 


Table  I - Dispersion  relation  polynomial  coefficients 

2 3 

f = a„  + a,  k + a-,k  + a-k 
0 1 2 3 


Specimen 

^0 

, 

(X  10-1) 

^2 

(X  10*3) 

^3 

(X10“°) 

i = 1 

(0.335  cm) 

0.2819 

0.1301 

0.3632 

-0.6759  i 

< 

i = 2 

(0.543  cm) 

0.2084 

0.1527 

0.3522 

-0.6160  j 

i = 3 
(0.878  cm) 

0.2567 

0.1270 

0.3788 

-0.7136  j 

i = 4 

(1.753  cm) 

0.2115 

0.1530 

0.3438 

j 

] 

-0.5604 

i = 5 

(2.569  cm) 

0.1935 

0.1569 

0.3362 

-0.5509 

i 

i 


A - 16  - 


I 


I 


Fiqure  Captions 


Figure  1 - Schematic  of  the  experimental  setup. 


Figure  2 - Measured  transducer  and  electronics  phase  delay.  Method  I; 
Transducers  in  contact;  Method  II:  Measured  using  a 
plexiglass  sample. 


Figure  3.  Original  phase  delay  data  obtained  from  shear  waves  in 
Boron-Epoxy,  propagating  along  the  fiber  direction,  with 
polarization  normal  to  the  ply  layers.  N*  is  the  cycles 
of  phase  relative  to  any  first  in-phase,  out-of-phase 
condition. 


Figure  4 - Phase  shift  of  the  specimens  i = 2, 3, 4, 5 for  various  phase 
shift  values  of  specimen  i = 1,  based  on  the  data  shown  in 
Figure  3. 

Figure  5 - Summed  phase  shift  obtained  from  Fig.  4 for  various  values  of 
the  phase  shift  for  the  data  of  specimen  i = 1. 

^ Figure  6 - The  dispersion  relation  derived  from  the  data  shown  in 

Figure  3. 

Figure  7 - The  corrected  dispersion  relation  data  for  the  0.543  cm  thick 
specimen.  Also  shown  is  the  polynomial  curve  fit  and  the 
derived  phase  and  group  velocity.  Open  circles  represent  the 
j group  velocity  determined  with  the  four-point  Lagrangian 

interpolation  technique. 

Figure  8 - Phase  and  group  velocities  obtained  from  all  the  data  shown  in 
Figure  3.  Comparison  is  made  with  the  results  obtained  from 
r.f.  burst  measurements  by  Tauchert  and  Guzelsu  (Ref.  11). 
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ON  THE  DETERMINATION  OF  PHASE  AND  GROUP  VELOCITIES 
OF  DISPERSIVE  WAVES  IN  SOLIDS 

Wolfgang  Sachse  and  Yih-Hsing  Pao 
Department  of  Theoretical  and  Applied  Mechanics 
Cornell  University,  Ithaca,  N.Y.  14853 

ABSTRACT 

A new  technique  is  developed  to  determine  the  dispersion  relation 
and  the  propagational  speeds  of  waves  in  dispersive  solids.  The  phase 
spectrum  of  a broadband  pulse  is  linearly  related  to  the  dispersion 
relation  of  the  dispersive  medium.  The  method  is  simpler  than  the  con- 
tinuous wave-phase  comparison  technique.  Application  is  made  to  mea- 
sure the  phase  and  group  velocities  of  waves  in  fiber-reinforced  com- 
posite materials  and  in  thin  wires.  This  technique  is  expected  to  be 
applicable  to  measurements  of  acoustic  or  electromagnetic  wave  speeds  in 
other  dispersive  media. 


I.  INTRODUCTION 


The  velocity  of  a wave  is  usually  associated  with  the  phase  difference 
between  the  vibrations  observed  at  two  different  points  during  the  passage 
of  the  wave.  A plane  harmonic  wave  of  amplitude  A,  angular  frequency  ui 
(=  2'irf)  and  wavelength  X(=  2TT/k),  propagating  in  a non-attenuating  medium 
may  be  represented  by  A cos((jjt  - kx+C)  = A cos(j(t  - kx/oj+e)  where 
is  a constant.  The  phase  of  the  wave  is  then  defined  by  the  argument  of  the 
cosine  function  and  the  wave  velocity  is  v = w/k,  which  is  usually  called 
the  phase  velocity. 

When  the  phase  velocity  v of  a plane  harmonic  wave  is  a function  of 
its  frequency  or  wavelength,  then  the  propagating  medium,  according  to 
Havelock  [1],  is  dispersive.  Such  dispersion  may  be  caused  by:  (i)  the 
presence  of  specimen  boundaries  (geometric  dispersion),  (ii)  the  frequency 
dependence  of  material  constants,  such  as  mass  density,  elastic  moduli, 
dielectric  constants,  etc.  (material  dispersion),  (iii)  the  scattering  of 
waves  by  densely  distributed,  fine  inhomogeneities  in  a material  (scattering 
dispersion),  (iv)  the  absorption  or  dissipation  of  wave  energy  into  heat,  or 
other  forms  of  energy  in  an  irreversible  process  (dissipative  dispersion),  and 
(v)  the  dependence  of  the  wave  speed  on  the  wave  amplitude  (non-linear 
dispersion).  We  shall  omit  the  last-mentioned  source  of  dispersion  in  the 
analysis  to  be  presented.  The  theory  of  the  propagation  of  pulses  in 
dispersive  media  is  well  known  [1-3].  The  principal  feature  is  that  the 
pulse  does  not  retain  its  initial  shape  as  it  propagates  through  the  dispersive 
medium.  Thus,  a short  duration  pulse  may  be  dispersed  into  wave  trains  in 
time. 
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Measurements  of  the  phase  versus  frequency  of  continuous  waves  in 
dispersive  media  determines  the  group  velocity  of  the  wave.  By  definition, 
the  group  velocity  U equals  du/dk  = v + k dv/dk.  From  the  data  of  U 
versus  u,  it  seems  that  one  could  evaluate  v as  a function  of  k by 
integrating  the  preceeding  equation.  However,  in  order  to  complete  the 
integration,  it  would  require  precise  values  of  u and  k at  the  lower 
limit  of  integrations  and  these  could  not  be  easily  obtained  experimentally. 

For  unidirectional  waves  in  a solid,  the  dispersion  characteristics 
are  similar  to  those  of  the  propagation  of  electromagnetic  signals  in 
transmission  lines,  or  those  of  ultrasonic  signals  in  delay  lines.  The  former 

is  disscu.sed  in  many  standard  texts  [4]  and  tho  latter  was  addressed  by  Young  [5]. 
Young  clei'^ly  defined  the  concept  of  "phase  delay"  of  continuous  waves,  and 

"group  dei.?y"  of  a narrow  bandwidth  pulse,  and  he  showed  that  the  latter 
is  related  lO  the  derivatives  of  the  phase  characteristics  of  the  pulse 
with  respe;i  to  frequency.  From  the  "group  delay,"  the  group  velocity  of 
a pulse  can  be  determined.  However,  Young  [5]  considered  only  waves  without 
attenuation,  and  pulses  with  narrow  bandwidth. 

Previous  experimental  work  in  solids  has  dealt  with  both  continuous- 
waves  and  pulses.  As  mentioned  above,  the  phase  velocity  may  be  frequency- 
dependent.  Thus  errors  in  its  determination  are  minimized  by  using  mono- 
chromatic waves.  In  the  so-called  7r-point  phase  comparison  technique,  which 
appears  to  have  been  first  used  by  Balamuth  [6]  and  is  recently 
reviewed  by  Papadakis  [7],  the  half-wavelength  of  the  wave  is  determined 
directly  by  varying  the  propagation  distance  between  source  and  receiver  and 
observing  the  corresponding  TT-phase  change  of  the  received  signal.  Although 
such  measurements  are  easily  carried  out  in  gaseous  and  fluid  media,  this  is 
not  the  case  in  solids. 
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The  main  difficulty  of  the  -rr-point  phase  technique  is  that  the  phase 
must  be  measured  to  a frequency  low  enough  so  that  the  phases  between  the 
transmitted  and  received  signals  differ  exactly  by  it,  or  a known  multiple 
of  IT  [7].  To  overcome  this  difficulty,  Lynnworth  et  al . [8]  used  a sliding 
wedge  technique  to  vary  the  propagation  length  in  3-dimensional  carbon- 
phenolic  composite  specimens  to  determine  the  wavelength  and  hence  the  phase 
velocity  directly.  Unfortunately,  the  wedge  interface  and  the  pronounced 
material  microstructure  can  give  rise  to  anomalous  effects.  Other  researchers 
used  specimens  of  fixed  thickness  and  varied  the  frequency  of  excitation. 

The  group  velocity  is  thus  determined  as  a function  of  frequency,  as  will 
be  discussed  in  nv^re  detail  in  Section  II,  Recently,  however.  Ting  and 
Sachse  [9]  have  combined  the  variable  frequency  method  and  the  variable  path- 
length  method.  They  developed  a constraint  equation  from  which  the  number  of 
wavelengths  in  a specimen  can  be  established  unambiguously. 

Pulse  measurements  of  the  dispersion  have  been  made  with  broadband 
pulses  as  well  as  narrowband  bursts.  In  a broadband  pulse  technique,  the 
frequency  dependent  reflection  coefficient  of  a pulse  from  a buffer-specimen 
interface  is  measured,  from  this,  the  frequency-dependent  phase  velocity 
of  the  specimen  inferred.  The  technique  has  been  applied  to  phase  velocity 
measurements  in  teflon  [8j  and  boron-epoxy  specimens  [10].  There  is  serious 
question,  however,  whether  the  measured  phase  velocity  is  that  which  charac- 
terizes the  bulk  material,  rather  than  the  interfacial  layer  at  the  buffer 
and  specimen  interface. 

In  burst  measurements,  one  measures  the  delay  of  a burst  with  a 
certain  center  frequency  through  the  thickness  of  a specimen.  Tauchert  [11], 
utilizing  pairs  of  transducers  with  differing  center  frequencies,  measured  the 
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group  velocity  in  this  way  in  boron-epcxy  composite  materials.  Using  a single 
pair  of  broadband  transducers  and  varying  the  frequency  of  the  electrical 
excitation  burst,  Sutherland  and  Lingle  [12]  measured  dispersive  effects  in 
boron-aluminum  specimens.  VJith  burst 

delay  measurements  in  very  dispersive  materials,  it  is  possible  to 
unambiguously  identify  equivalent  points  in  the  undelayed  burst  and  the 
specimen-delayed  burst  therefore,  it  is  not  clear  what  wavespeed  is  actu- 
ally measured  [7].  Furthermore,  as  emphasized  by  Young  [5],  the  envelope 
of  a received  signal  may  differ  considerably  from  the  input  signal  and  this 
makes  it  difficult  to  determine  its  delay-time  through  a specimen.  This 
problem  can  be  minimized  by  amplitude  modulating  the  r.f.  carrier  with  a 
particular,  smooth  envelope  such  as  a Gaussian  contour  [5], 

Other  examples  of  elastic  wave  dispersion  measurements  are  those  with 
a shock  tube  technique  by  Whittier  and  Peck  [13]  ,in  carbon-phenolic  composites; 
by  Yew  and  Yogi  [14]>in  steel-epoxy  specimens  with  an  optical  interferometry 
technique;  and  by  Kaarsberg  [15]  in  geological  materials  with  ultrasonic 
techniques. 

Our  investigation  was  motivated  primarily  by  the  determination  of  wave 
speeds  in  eutectic  and  fiber-reinforced  composite  materials,  and  porous  solids. 
In  Section  III  of  this  paper  we  review  the  principles  of  the  ir-point  phase 
comparison  method  and  its  limitations  in  uniquely  determining  the  dispersion 
relation  of  a material. 

In  Section  IV,  we  propose  a different  technique  - based  on  the  phase 
spectral  analysis  of  a broadband  pulse  - for  the  determination  of  phase  and 
group  velocities  of  wa/es  in  dispersive  media.  In  this  technique,  the 
phase  function  of  a Fourier-analyzed  pulse  is  evaluated.  It  is  shown  how 
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the  phase  and  group  velocities  as  well  as  the  dispersion  relation  for  the 
material  can  be  determined  directly  from  it.  The  experimental  technique 
and  its  implementation  are  described  in  Section  V.  In  Section  VI  we  show 
the  results  obtained  with  this  technique  to  characterize  dispersive  and 
non-dispersive  materials  and  to  compare  them  with  the  results  obtained  with 
the  continuous  wave,  ir-point  phase  comparison  technique. 

Although  the  examples  we  cite  in  this  paper  are  limited  to  ultrasonic 
stress  waves  and  pulses  propagating  in  dispersive  solids,  the  principle  and 
the  experimental  technique  described  in  Sections  III  and  IV  are  perfectly 
general,  such  that  similar  observations  can  be  made  with  linear  acoustic, 
electro-magnetic,  or  optical  waves  propagating  in  dispersive  media.  A 
general  discussion  of  the  dispersion  of  waves  and  the  concept  of  group 
velocity  can  be  found  in  Refs.  [16]  and  [17]. 

II.  DISPERSION  OF  A LINEAR,  CAUSAL  SYSTEM 

A solid  specimen  sustaining  the  propagation  of  a plane  wave  along  a 
length  or  thickness  dimension,  comprises  a one-dimensional,  linear  system. 
Such  systems  admit  plane  wave  solutions  of  the  form: 

u(x.t)  = (1) 

As  noteo  earlier,  the  phase  velocity  v(=  oj/k)  is  a constant  for  non- 
dispersive  media.  More  generally,  v is  a function  of  w.  Also,  in 
real  systems,  the  wave  amplitude  gradually  attenuates  through  beam 
spreading,  scattering,  mode  conversion,  coupling  losses  between  transducer 
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and  specimen  and  the  absorption  and  dissipation  of  energy  in  the  specimen 
material.  As  long  as  the  total  energy  is  conserved,  w must  be  real  and 
the  attenuation  of  the  wave  amplitude  can  be  accounted  for  by  letting  k 
be  complex.  That  is 

k = B + ia  (2) 

For  waves  propagating  in  the  direction  of  +x,  Eq.  (1)  is  then  written  as 

u(x,t)  = (3) 

In  general,  both  a and  B in  Eqs.  (2)  and  (3)  are  functions  of  u. 

Although  the  function  in  (3)  is  not  periodic,  we  can  still  define  a 
phase  velocity  v,  and  group  velocity  U of  the  wave  as 


V 


i»  - w.  I - do) 
Reiry  ■ B ’ dB 


(4) 


The  wave  length  of  the  decaying  sinusoidal  wave  is  2tt/B,  and  a is  the 
attenuation  factors.  We  caution  that  the  group  velocity  U as  customarily 
defined  only  has  physical  significance  when  a(a))  is  small  [16].  By  small 
a(a))  we  mean  that  the ‘attenuation  in  nepers  for  one  wave  length  is  small. 

The  dispersive  property  of  a physical  system  is  characterized  by 
specifying  the  a and  b as  functions  of  w.  For  a causal  system,  in 
which  the  response  of  the  system  to  a driving  force  cannot  precede  the 
first  arrival  of  the  signal  generated  by  the  driving  force,  the  aiu)  and 
6(a))  are  related  by  Kramers-Kronig  relation.  Mathematically,  for  a causal 
system,  the  ,j(a))  is  related  to  the  Hilbert  transform  of  B(a))»  and 
vice  versa  [3].  Thus  we  shall  regard  that  aU)  can  be  calculated  from 
B(a)),  and  the  dispersion  property  of  the  linear  causal  system  is  defined 
by  the  relation 
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B = 6(u),  or  CO  = u)(B)  (5) 

^ Ak 

In  Eq.  (5),  B((o)  is  a real  function  of  the  real  variable  co,  and  aj(B) 

is  the  inverse  function  of  B(co).  The  objective  of  the  experiments  is  to 
determine  the  dispersion  relation,  Eq.  (5),  for  a solid  medium,  and  then  to 
calculate  v((o)  and  U((o)  from  Eq.  (4);  or  to  measure  either  v or  U and 
then  determine  the  dispersion  relation,  Eq.  (5). 

If  there  exists  a dissipative  mechanism  in  the  system,  the  frequency 
(0  may  also  be  complex  valued.  One  example  is  the  thermoelastic  wave  in 
solids  [18,  19].  In  such  a case,  there  are  two  choices:  one  is  to 
restrict  the  answer  to  real  co  but  complex  k,  the  other  is  to  limit  the 
answer  to  real  k but  complex  co.  The  corresponding  phase  velocities  are 
defined  respectively  as  v = aj/Re(k)  and  v = Re(oj)/k,  and  the  two  results 
do  not  always  agree.  The  meaning  of  group  velocity  also  becomes  uncertain. 
The  situation  is  further  complicated  by  possible  anisoti'opy  in  the  material. 
The  concept  of  phase  and  group  velocity  for  waves  in  a dissipative  and 
anisotropic  medium  has  been  discussed  in  Ref.  [19], 

III,  THE  METHOD  OF  m-POINT  PHASE  COMPARISON 

By  transmitting  continuous  waves  of  varying  frequencies  through  a 
specimen  of  length  L,  a phase  versus  frequency  relation  can  be  measured 
for  the  specimen.  By  adjusting  the  frequency  to,  say,  f^,  the  signal  at 
the  receiver  is  first  brought  into  in-phase  or  m-radians  out-of-phase 
with  that  from  the  transmitter.  At  this  frequency,  let  the  wavelength 
be  X.J,  and  the  number  of  complete  sine  waves  by  thus  f^  = where  v^ 

is  the  phase  velocity.  An  increment  of  the  phase  by  a known  multiple  of 
n radians  is  made  by  changing  the  frequency  to  with  number  of  sine 
waves  of  length  a^,  and  f2  = VgNg/L. 
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For  nondispersive  waves,  = V2  = v.  By  setting  N2  = + 1 in 

the  experiment  [7],  we  find 


V = L(f2-f^) 


(6) 


In  this  situation,  it  is  not  necessary  to  know  the  exact  value  of  at 
the  starting  point. 

For  dispersive  media,  v^  differs  from  because  of  the  change  of 
frequency,  and  so  one  cannot  determine  either,  unless  is  known  exactly. 
However,  as  illustrated  by  the  hypothetical  experimental  data  shown  in  Fig. 

1 , it  is  still  possible  to  determine  the  group  velocity  U by  measureing 
= 1,2,3,...)  and  N.  where  = N.  + 1.  This  follows  from  the 

following  relations; 


U = (L  ^)  = L 


(7) 


In  this  case,  U is  only  determined  in  a range  of  frequencies  higher  than  f^ . 

Once  U(N)  is  determined  for  a range  of  ^ , the  phase  velocity 

can  be  calculated  by  the  following  integration 

1 ^ 

V = ^ ; U(N)dN  + V,  — (8) 

' H 

Thus,  we  are  facing  an  uncertainty  unless  the  value  for  is  known.  This 
can  easily  be  seen  from  Fig.  1,  in  which  the  location  of  is  uncertain.  The 
slope  of  the  tangent  to  the  curve  f vs.  N which  is  related  to  the  group  velocity 
can  be  calculated  for  N > but  the  slope  of  the  secant  (a  line  from  the  origin 
to  a point  on  the  curv?’)  cannot  be  determined.  The  slope  of  the  secant  is 
precisely  v/L. 
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IV.  THE  METHOD  OF  PHASE  SPECTRAL  ANALYSIS 
Fourier  Synthesis  of  the  Transmitted  Pulse 

It  was  mentioned  that  a linear  system  admits  plane  harmonic  waves  des- 
cribed by  Eq.  (1),  and  if  attenuation  exists,  by  Eq.  (3).  Based  on  the 
principle  of  superposition,  a pulse  u(x,t)  propagating  in  the  linear 
medium  x > 0 may  be  expressed  as  an  addition  of  all  plane  harmonic  waves. 


u(x,t)  = ^ / do)  / 


(9) 


where  A(c)  is  an  unspecified  amplitude -function,  and  both  v and  a are  functions  of 
Assuming  that  the  medium  is  at  rest  initially,  and  a disturbance  F(t)  is 
generated  at  the  location  x = 0,  for  t ^ 0,  that  is 

u(0,t)  = F(t)  (10) 

Then,  from  Eq,  ( 9),  we  find 


F(t)  = ^ 7 do,  7 A(Oe'"^^-^^dc 


(11) 


This  expression  is  exactly  the  Fourier  integral  representation  of  the  input  time 
function  F(t).  Hence  A(c)  = F(5),  ana  the  amplitudes  of  the  superposed  plane 
waves  are  fixed  by  the  end  condition  at  x = 0.  We  thus  write  Eq.  (9)  as 

u(x,t)  = ^ f <L  f F(c)e’“’‘e’“^^’’^'''’''^^dc  . (12) 

• OJ  *00 

If  a and  V are  constant,  we  find  from  Eq.  (12)  and  the  theorem  of 
Fourier  integral  the  well-known  result 

u(x,t)  = F(t-x/v)e‘'‘'‘  (13) 

The  input  pulse  F(t)  at  x » 0 propagates  through  the  non-dispersi ve  medium 
without  a change  of  shape,  except  for  an  exponential  decrease  in  amplitude. 
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When  both  a and  v are  functions  of  cd,  the  double  integral  in 
Eq  . (12.)  gives  rise  to  the  dispersed  pulse  at  any  location.  To  determine 
the  signal  at  x = L in  a specimen,  we  evaluate  first  the  integration  with 
respect  to  C, 


where 


u(L.t) 


F(aj)  = 


du 


(11*) 


F(c)e"''‘^dC  (15) 

Note  that  F(co)  is  the  Fourier  transform  of  the  input  F(t)  which  is 
supposed  to  be  known.  Since  B = w/v,  Eq.  (l4)  can  be  rewritten  as 

r 00 


u(L,t)  = 2tt 


r^/  \Q”^-"i6L-i  iujt. 

LF(u))e  e je  dw 


(16) 


The  exact  shape  u(L,t)  can  be  calculated  from  the  above  integral  if 
the  dispersion  relation  (Eq  • 5 ) is  known,  but  the  calculation  will  be  in 

general  be  very  difficult.  However,  the  spectrum  of  the  transmitted  pulse  at 
X = L is  given  by  the  Fourier  transform  of  u(L,t),  which  is  precisely 
the  quantity  inside  the  brackets  of  Eq  . (i6). 


u(L,u)  = F(u))e  *^6 


(IT) 


In  general,  the  Fourier  transform  of  a causal  time  function  F(t)  may  be 
complex.  Thus 

F(a,)  = |F(cu)|e"''‘^  = |F(a,)le  (l8) 


where  4)^  = ojC*  Eqi.  (17)  can  then  be  written  as 


17(1,0))  = |F(o))e‘“‘-|e‘‘'^®'-'’‘^^ 


(19) 


This  means  that  the  phase  spectrum  of  u(L,t)  is  given  by 

- — - — • 


- B r^ll 


A 

(J>(u)  = BL  + (i)^  = L8(ti))  + ojC 


It  is  seen  from  Eq  . (20)  that  the  phase  spectrum  of  the  transmitted  pulse  is 
linearly  related  to  the  dispersion  relation  6(01)  of  the  meidum. 

Phase  Spectra  and  Dispersion  Relations 

In  the  next  section,  we  show  experimentally  that  an  accurate  spectrum 
of  the  transmitted  pulse  u(L,t)  can  be  determined  by  digital  analysis. 

The  spectrum  is  composed  of  two  parts:  the  amplitude  spectrun  and  the  phase 
spectrum.  The  latter  is  precisely  the  (})(u))  of  Eq.  (20).  Thus,  the 
dispersion  relation  of  a medium  is  determined  uniquely  by  the  phase  spectrum 
of  a propagating  pulse, 

3(u))  = ^ C(1)(u>)-<)1q]  (21) 


The  phase  velocity  is  obtainable  from  the  experimental 


data. 


v(a))  = ^ 


ojL 


3 


(22) 


By  differentiating  the  phase  function,  we  obtain 

A 

= I ^ + £ = - L _ + £ = k + r 

do)  dcj  ^ dcj/d6  U 


(23) 


The  group  velocity  U can  then  be  calculated  from  the  derivatives  of  the 
phase  spectrum 


u ■ L(^  - 5) 


(24) 


When  (D  is  complex-valued,  as  in  a dissipative  medium,  the  Fourier 
synthesis  in  Eq.  ( 5)  is  meaningless  unless  the  path  of  integration  for  the 
complex  variable  ui  is  specified.  An  alternative  would  be  to  assume  k 
real  valued  and  to  write 


u(x,t)  = 


1 

iSL 


dk 


(25) 
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The  c is  a phase  constant  in  the  space-axis,  and  v = w/k  must  be 
complex  valued  for  a medium  with  damping.  As  shown  by  Havelock  [1],  an 
integral  representation  such  as  Eq.  (25)  is  suitable  for  analyzing  initial 
value  problems  (u  is  given  at  t = 0).  For  waves  in  a medium  without 
dissipation,  these  two  representations  (Eqs,  9 and  25  ) are  equivalent. 


A Simple  Example 


That  the  phase  spectrum  of  a propagating  pulse  is  related  to  the  group 
velocity  may  be  illustrated  by  the  following  simple  example.  Consider  a 
square  pulse  of  duration  21^  generated  at  x = 0, 

u(0,t)  = Uq  0 < t < 2Tg  (26) 


and  u(0,t)  =0  for  t < 0,  and  for  t > 21^.  At  time  t = T,  the  pulse 
travels  through  a non-dispersive  and  non-attenuating  medium  a distance  L. 
Since  the  pulse  shape  is  preserved  (Eq.  13),  the  disturbance  at  x = L is 

u(L.t)  = Uq  T < t < T+2Tq  (27) 


and  u(L,t)  * 0 for  t < T and  t > T+2Tjj.  The  Fourier  transform  of  u(0,t) 
and  u(L,t)  are  respectively 


u(.0,aj)  = 2Tq 
17(1,0))  = 2Tq 

The  phase  spectrum  of  7(1,0)) 


sino)T^  “iwT^ 
0 „ 0 


sinojT^  -io)(T+TQ) 
and  its  derivative  are 


(28) 

(29) 


4)(cj)  = a}(T+T^),  d())/du)  = T+T^ 


i 

\ 


i 


I 


(30) 
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Comparing  Eq.  (29)  with  Eq.  (19)  with  a = 0,  Eq.  (30)  with  Eq.  (23). 
we  identify  C = Tq  and  find  T = L/U.  This  last  result  states  that 
the  square  pulse  traverses  the  distance  L with  the  group  velocity  U 
which,  in  this  case,  is  also  equal  to  the  phase  velocity. 

V.  EXPERIMEns  AND  DIGITAL  SPECTRAL  ANALYSIS 

The  two  methods  described  in  the  previous  sections  were  implemented  to  measure 
the  dispersion  relation  and  phase  and  group  velocities  for  non-dispersive  and 
dispei^sive  materials.  The  experimental  set-up  for  the  continuous  wave 
iT-point  phase  difference  measurements  is  similar  to  that  in  Ref.  7 (Fig. 

22).  The  block  diagram  for  the  equipment  used  in  the  pulse  transmission 
experiments  is  shown  in  Fig.  2.  Two  nearly  identical  broadband  • 
transducers  were  used  in  both  experiments. 

In  the  continuous  wave  experiments,  the  source  transducer  was  excited 
with  a 10  volt  peak-to-peak  harmonic  wave  ranging  in  frequency  from  50  kHz 
to  20  MHz.  The  signals  from  the  receiver  transducer  were  amplified  and 
displayed  on  an  oscilloscope.  The' ampl itude  and  phase  measurements  of  the 
continuous  waves  were  made  directly  from  the  oscilloscope,  and  with  a 
vector  voltmeter  for  frequencies  above  1 MHz. 

In  the  pulse  experiments,  the  source  transducer  was  excited  with  a 
-250  volt  pulse  of  35  nsec  duration.  A sampling  oscilloscope  was  used  to 
overcome  the  bandwidth  limitation  of  the  A/D  converter  attached  to  a mini- 
computer. Any  selected  portion  of  a time  record  could  be  sampled  by  means 
of  a digital  delay  unit  used  in  conjunction  with  the  sampling  plug-in  unit 
of  the  oscilloscope. 
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A Fast-Fourier  transform  algorithm  was  used  to  compute  the  Fourier 

amplitudes  of  the  input  time  function  F(t).  The  actual  time  record  is 

associated  with  1024  time  domain  data  points  which  were  padded  with  1024 

zeroes.  The  real  and  imaginary  parts  of  the  Fourier  transform,  F^(u)) 

and  ^2(0)).  respetively,  were  first  calculated.  The  absolute  value  of 

-1/2 

the  transform  is  given  by  jF(uj)l  = [F-j  (ui)F2(id)  j , and  the  phase  func- 
tion by  = tan'^[F2(uD)/F^(w)]. 

The  computer  algorithm  limits  the  values  of  $ in  the  reduced  range  from 
TT  to  -IT.  Thus,  discontinuities  of  it  for  the  phase  spectra  were  found 

whenever  the  magnitude  of  (})(a))  becomes  less  than  or  greater  than  zero  and 
2n.  This  was  corrected  for  by  adding  or  subtracting  the  amount  of  jump  to 
render  the  phase  spectrum  continuous.  The  resultant  data  was  then  used  to 
compute  the  wavenumber,  phase  velocity  and  group  velocity  associated  with 
the  pulses  u(0,t)  and  u(L,t)  according  to  Eqs,  (21),  (22)  and  (24). 

To  initialize  the  pulse  analysis  system  and  thus  set  <{)^,  the  source 
transducer  was  placed  in  intimate  contact  with  the  receiving  transducer  and 
the  maximum  of  the  received  excitation  pulse  was  made  to  coincide  with  the 
start  of  the  sampled  time  record  by  means  of  the  vernier-controlled  digital 
delay  unit.  In  this  way  was  set  equal  to  0.  Alternatively,  the  phase 
spectrum  of  the  received  excitation  pulse  could  be  computed  and  stored  in 
the  data  analysis  system.  This  spectrum  is  exactly  and  it  could  be 
subtracted  from  all  subsequent  phase  spectra  calculations. 
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The  sampled  sweep  was  digitally  delayed  in  steps  of  100  nsec  until 

the  pulse  to  be  analyzed  was  displayed  on  the  sampling  oscilloscope.  The 

time  shift,  x , in  increments  of  100  nsec,  relative  to  x = 0.0  was 
s s 

read  into  the  computer  as  was  the  propagation  distance,  L.  Equations 
(21 ).  (22)  and  (24)  for  the  dispersion  relation,  phase  and  group  velocity, 
respectively  were  modified  to  include  the  time  shift  factor.  They  are, 
in  terms  of  frequency  f, 

^f)  - (|)(f)]  + 27rfXg  (51) 


v(f)  = 


!4>q  - + 27rfXg 


U(f)  = 27rL/(d4,/df-Tg) 


The  operation  of  the  system  was  checked  by  entering  a simulated  square  pulse, 
and  the  phase  spectrum  ())((d)  was  indeed  a straight  line  as  given  by  Eq.  (30). 


VI.  EXPERIMENTAL  RESULTS 


To  experiment  with  waves  in  a non-dispersive  medium,  a broadband 
longitudinal  ultrasonic  pulse,  comprised  mainly  of  frequency  components 
from  3 to  12  MHz  was  propagated  through  a 1.900  cm  thick  specimen  of  6061 -T6 
aluminum.  In  Figure  3(a)  is  shown  the  excitation  pulse  u(0,t)  and  in 
Fig.  3(b)  the  first  received  pulse  u(L,t)  at  L = 1.900  cm.  In  Figures 
3(c)  and  3(d)  are  shown  the  Fourier  amplitude  spectrum  and  phase  spectrum 
of  u(L,t).  Fig.  3(d)  also  yields  the  dispersion  relation  as  B = (})/L 
(Eq.  21  ).  The  calculated  phase  velocity  and  group  velocity  of  this  specimen 
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are  shown  in  Figs.  3(e)  and  3(f),  It  is  clear  from  Fig.  3(d)  that  there 
is  little  dispersion  in  this  material  in  the  frequenty  range  from  0 to  20 
MHz. 

The  circular  dots  shown  in  Figure  3(f)  represent  the  group 
velocity  values  determined  from  the  data  ij';easured  by 
the  continuous-wave,  n-phase  technique  (Eq.  7).  For  this  specimen,  no 
measurements  were  possible  below  1.5  MHz  because  transducer  near-field 
effects  dominated  the  observations.  At  frequencies  above  7 MHz,  data 
were  recorded  only  at  integer  frequency  points.  The  average  group  velocity 
measured  using  the  continuous-wave  technique  was  0.618  + 0.009  cm/psec. 

The  results  for  a dispersive  medium  such  as  composite  materials  are 
shown  in  Figure  4.  A broadband  shear  pulse  is  propagated  through  a 0.546 
cm  thick  specimen  of  96-ply  Boron-Epoxy.  In  this  case,  the  wave  propagates 
in  the  direction  of  the  fiber  and  the  particle  displacement 
direction  is  perpendicular  to  the  fiber  and  parallel  to  the  ply  layers. 

Figs.  4(a)  and  4(b)  show  the  initial  pulse  and  the  dispersed  pulse,  respectively, 
detected  at  the  receiving  transducer.  It  is  apparent  that  the  high  frequency 
components  of  the  pulse  propagate  faster  than  the  low  frequency  components. 

The  relative  amplitudes  of  the  various  frequency  components  in  the  pulse 
are  shown  in  Fig.  4(c).  As  expected,  the  high-frequency  components  are 
markedly  lower  in  amplitude  than  the  low-frequency  components.  The  dispersion 
relation  is  shown  in  Fig.  4(d),  which  is  no  longer  a linear  function  of  oj. 

This  curve  is  obtained  from  the  experimental  phase  spectrum  by  scaling  the 
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ordinate  with  the  length  L.  Figures  4(e)  and  4(f)  show  the  results 
of  the  phase  and  group  velocity  computation  for  the  frequency  range  0 
10.0  MHz.  Results  for  the  spectral  amplitudes  beyond  10.0  MHz  are  not 
dependable,  as  the  amplitudes  become  negligibly  small  there. 

The  circular  dots  on  the  group  velocity  curve  are  obtaineo  from  the 
ir-phase  technique.  The  extent  of  the  dispersion  can  be  readily 

seen  in  the  dispersion  curve  in  Fig.  4(d). 

As  another  demonstration  of  the  application  of  the  Fourier  phase 
analysis  technique,  we  show  its  use  in  the  analysis  of  a broadband  pulse 
propagating  in  a 76.4  cm  length  of  1.61  mm  (OD)  Remendur  tubing.  In  this 
case,  both  source  and  receiving  transducers  consisted  of  280  turn  coils  of 
#32  AWG  wire  wound  on  a section  of  glass  tubing.  The  coils  were  slipped 
over  the  ends  of  the  Remendur  tubing  and  a -50  volt  pulse  having  a width 
of  about  1 ysec  was  used  as  excitation.  The  incident  pulse  contains 
frequency  components  from  150  to  500  KHz.  The  signal, received  with  a 
second  coil  placed  over  the  other  end  of  the  tubing,  is  shown  after  40  db 
amplification  in  Fig.  5(a).  The  pulse  has  been  shifted  by  135.7  ysec  to 
bring  it  on  scale.  In  contrast  to  the  previous  example,  the  low  frequency 
components  of  the  pulse  in  this  case  propagate  faster  than  the  high 
frequency  components.  The  amplitude  spectrum  normalized  to  the  maximum 
value  is  shown  in  Fig.  5(b),  while  in  Fig.  5(c)  are  shown  the  derived 
dispersion  relation,  phase  and  group  velocities  for  this  material. 


VII.  CONCLUSIONS 

In  this  paper,  it  has  been  shown  that  the  pulse  phase  spectrum  de- 
fines the  dispersion  relation  of  a linear  causal  medium,  ana  a phase 
spectral  analysis  of  a dispersed  pulse  yields  results  for  both  phase  and 
group  velocities  of  waves  in  dispersive  media. 

Most  continuous-wave  techniques  require  tedious  frequency  scanning 
and  oata  reauction.  In  contrast,  the  pulse  analysis  technique  is  rapid. 

The  technique  appears  not  to  be  restricted  to  the  analysis  of  elastic 
pulses  propagating  in  solids.  It  could  well  find  application  to  studies  of 
the  dispersion  of  acoustic  or  electro-magnetic  waves  in  dispersive  media. 
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Figure  Captions 


Figure  1 
Figure  2 
Figure  3 


Figure  4 


Figure  5 


Hypothetical  dispersion  curve-frequency  (f)  versus  number 
of  waves  in  a length  L(N). 

Block  diagram  of  the  apparatus  for  measuring  the  phase 
spectrum. 

Pulse  propagation,  dispersion,  phase  and  group  velocities  of 
waves  in  6061 -T6  Aluminum 

(a)  Input  P-  pulse  in  voltage  (tj  = 0.0  usec). 

(b)  Received  P-pulse  at  L in  voltage  (t^  = 2.6  usec). 

(c)  The  amplitude  spectrum  of  the  received  pulse  (Relative 
amplitude  versus  frequency  f in  mega  hertz). 

(d)  The  phase  spectrum  of  the  received  pulse  and  the  disper- 
sion relation  (wave  number  B in  1/cm  versus  frequency 
f in  mega  hertz). 

(e)  Phase  velocity  v (cm/ysec)  versus  frequency  f 
(mega  hertz). 

(f)  Group  velocity  U (cm/ysec)  versus  frequency  f 
(mega  hertz). 

Pulse  propagation,  dispersion,  phase  and  group  velocities  of 
waves  in  96-ply  boron-epoxy  composite  materials.  (Shear  wave 
propagating  in  the  direction  of  the  fibers) 

(a)  Input  S-pulse  in  voltage  (t^  = 0.0  ysec). 

(b)  Received  S-pulse  at  L in  voltage  (t^  - 0.0  ysec). 

(c)  The  amplitude  spectrum  of  the  received  pulse  (Relative 
amplitude  versus  frequency  f in  mega  hertz). 

(d)  The  phase  spectrum  of  the  received  pulse  and  the  disper- 
sion relation  (wave  number  6 in  1/cm  versus  frequency 
f in  mega  hertz). 

(e)  Phase  velocity  v (cm/ysec)  versus  frequency  f (mega 

hertz) . 

(f)  Group  velocity  U (cm/ysec)  versus  frequency  f (mega 

hertz) . 

Dispersion  of  a pulse  in  a Remendur  Tubing  of  1.61  mm  outer 
diameter 

(a)  Received  extensional  pulse  at  L in  voltage  (t^  = 135.7  ysec). 

(b)  The  amplitude  spectrum  of  the  received  pulse  (relative 
magnitude  versus  frequency  f in  mega  hertz). 

(c)  Dispersion  relation  (wave  number  B in  1/cm  versus  frequency 
in  mega  hertz),  phase  velocity  v in  cm/ysec,  and  group 
velocity  U in  cm/ysec. 
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ABSTRACT 

A scattering  matrix  approach,  that  involves  only  the  transition 
matrix  of  a single  obstacle,  is  proposed  for  studying  the  miiltiple 
scattering  of  scalar  waves  in  a medium  (matrix)  containing  identical, 
long,  parallel,  randomly  distributed  cylinders  of  arbitrary  cross 
section.  The  elastic  properties  of  the  cylinders  are  assumed  to  be 
different  from  those  of  the  matrix.  A statistical  approach  in  con- 
junction with  Lax's  "quasicrj'staJ.line"  approximation  is  employed  to 
obtain  equations  for  the  average  amplitudes  of  the  scattered  and 
exciting  fields  which  may  then  be  solved  to  yield  the  dispersion  rela- 
tions of  the  composite  medium.  Dynamic  elastic  properties  of  the 
composite  medium  containing  circular  and  elliptical  cylinders  are 
foxmid  in  the  Rayleigh  or  low  frequency  limit.  Nmerical  results  dis- 
playing phase  velocity  and  damping  effect  of  the  composite  medi\im 
are  presented  for  a wide  range  of  frequencies. 
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IMTRODUCTION 


The  subject  of  mviltiple  scattering  of  waves  is  of  interest  in 
many  fields  of  engineering  and  science.  In  acoustics,  it  has  important 
practical  applications  in  studies  of  distribution  of  flaws  in  solids, 
fiber  reinforced  composites,  porous  media,  fog  horn,  underwater  signal 
transmissions  etc.  Ever  since  Rayleigh's  analytical  treatment  of  scat- 
tering of  waves  by  randomly  distributed  particles  to  study  the  color 
of  the  sky^,  diffraction  sind  scattering  of  waves  have  been  the  subject 

of  extensive  research.  We  cite  here  papers  that  are  related  to  our 
2-11 

present  analysis  , in  which  additional  references  may  be  found.  Most 
of  the  available  results  are,  however,  confined  to  scatterers  of  simple 
geometries  like  circular  cylinder,  sphere  etc.,  and  often  at  low  frequencie 
A detailed  analysis  of  mviltiple  scattering  of  waves  by  obstacles  of 
arbitrary  shapes  and  over  a wide  range  of  frequencies  is  still  lacking. 

In  scattering  theories,  approximations  are  usually  made  at  a 
very  early  stage  for:  (l)  the  geometry  of  the  obstacle,  (2)  size  of 
the  obstacle  relative  to  the  wavelength  of  elastic  waves  and  (3'  dis- 
tribution of  the  obstacles  in  the  medi-um.  The  approximations  with 
respect  to  geometry  and  size  are  related.  If  the  obstacle  is  small 
compared  to  the  incident  wavelength,  it  is  not  possible  to  'see'  exact 
details  of  the  obstacle  and  usually  one  is  content  to  obtain  the  gross 
scattering  properties  of  the  obstacle.  This  is  the  so  called  Rayleigh 
limit  or  low  frequency  limit,  and  yields  corrections  to  the  solution 
for  point  scatterers.  As  far  as  the  distribution  of  the  obstacles 
is  concerned  one  either  has  regular  arrays  of  obstacles,  or  a random 
distribution.  In  the  latter  case  we  employ  a configurational  averaging 
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procecfure.  If  the  concentration  of  obstacles  is'  small,  i.e.,  they 
are  sparsely  distributed  we  may  employ  a single  scattering  or  first 
Born  approximation. 

Foldy^,  Lax^’^,  Twersky^’^,  Waterman  and  Truell”^  and  Karl  and 

g 

Keller  have  used  a statistical  approach  to  the  mxiltiple  scattering 

9 10 

of  scalar  waves.  Recently,  Eose  and  Mai  and  Mai  and  Chatterjee  have 
extended  the  statistical  approach  to  study  the  scattering  of  elastic 
waves  by  circular  cylinders  and  obtained  solutions  in  the  Rayleigh 
limit  and  under  static  loading,  respectively.  Datta^^  has  studied  the 
scattering  of  scalar  waves  by  cylinders  of  elliptical  cross  section 
using  the  method  of  matched  asymptotic  expansions  in  conjunction 
with  first  Born  approximation  in  the  Rayleigh  limit. 

The  objective  of  the  last  mentioned  papers  (Refe.  9>  10  and  11 ) 

is  to  determine  the  dj'namic  characteristics  of  a medivim  reinforced  by 

cylindrical  fibers.  The  same  objective  has  also  been  pursued  by  many 

authors  using  different  theories  and  a review  of  the  subject  is  given 
12  13  Ik 

by  Achenbach  , Peck  and  Moon  . Generally,  a dispersion  rela- 
tion for  plane  waves  in  the  composite  medium  is  determined  in  these 
papers,  but  seldom  is  a discussion  of  attenuation  or  damping  of  the 
medium.  In  the  Rayleigh  limit,  this  may  be  due  to  the  fact  that  the 
effective  wave  number  in  the  composite  medium  is  a real  quantity, 
and  in  the  case  of  periodic  layering,  the  periodic  nature  of  fibers  and 
the  boundedness  nf  the  body  in  a direction  normal  to  the  wave  vector 
confines  all  the  energy. 

On  the  other  hand,  Sutherland  and  Lingle^^  and  Tauchert  and 
Guzelsu^^  have  observed  appreciable  damping  in  these  materials  and  a 
few  experimental  values  have  been  given.  However,  they  also  stated  that 


■ the  actenuation  measurements  are  difficult  to  make  for  a vide  range 

of  frequencies. 


In  this  paper  ve  applj’’  a scattering  matrix  approach  to  multiple 

17 

scattering  of  horizontally  polarized  shear  vaves.  Waterman"^  intro- 
duced the  matrix  method  to  study  the  response  of  a single  obstacle 
to  acoustic  and  electromagnetic  vaves.  The  matrix  method  has  the 
advantage  of  being  applicable  to  obstacles  of  smooth  but  arbitrarj’-  shape 
and  also  at  vave-lengths  comparable  to  obstacle  size.  Peterson  and 
storm  have  given  formal  expressions  for  the  scattering  matrix  of  a 
distribution  of  obstacles  but  they  are  difficult  to  compute  in  practise. 


We  propose  an  alternate  method  vhich  combines  the  scattering  matrix  of  a 
I single  obstacle  and  a suitable  statistical  averaging  technique.  This 

approach  can  be  used  conveniently  to  obtain  explicit  solutions  at  the 
Rayleigh  limit,  and  to  obtain  the  dispersion  relation,  phase  velocity 
I and  attenuation  constants  of  inhomogeneous  media  for  a vide  range  of 

* frequencies. 

We  consider  a homogeneous,  isotropic  infinitely  extended  solid 
containing  K number  of  indentical,  long,  parallel,  randomly  distributed 
cylinders  of  arbitrary  cross  section.  We  assume  that  time  harmonic  plane 
shear  vaves  polarized  parallel  to  the  axis  of  the  cylinders  are  incident 
normal  to  the  scatterers  (SH  vaves).  When  H is  very  large,  a statis- 
tical analysis  is  presented  using  the  scattering  matrix  formulation. 

The  statistically  averaged  equations  are  then  solved  by  Lax's  "quasi- 
crystalline" approximation  to  yield  the  propagation  characteristics  of 
the  "average  waves"  in  the  medim.  Analj'tical  results  are  obtained 
for  the  wave  speed,  dispersion  relation  and  effective  elastic  moduli 
of  the  composite  in  the  Rayleigh  limit,  and  the5^  are  shown  to  agree  with 


the  existing  resvilts.  Numerical  results  are  also  presented  for  propagation 
speeds  and  attenuations  at  higher  frequencies. 


II.  TRANSITION  MATRIX  FORI-IULATION  FOR  I^^ULTIPLE  SCATTERING 

We  consider  N arbitrary  shaped,  long  and  parallel  cylindrical 

scatterers  with  piecewise  smooth  bovindary  surface  S embedded  in  an 

infinitely  extended  solid  (matrix)  which  are  referred  to  a coordinate  system 

as  shown  in  figure  la.  0^  and  0^  denote  the  centers  of  the  i"^^  and 

inclusions  and  they  can  be  represented  by  the  polar  coordinates  r^,  6^ 

and  z.  and  r,,  6.  and  z.,  respectively.  Let  p,  y be  the  density  and 
1 J j J 

rigidity  of  the  matrix  and  p^,  those  of  the  scatters. 

We  represent  an  incident  plane  SH  wave  of  unit  amplitude,  frequency 
0)  and  wave  vector  k by 


Q i(k*r-U)t) 

u=0,v  = 0,w  = w^  = e . (l) 

The  wave  is  assumed  to  move  along  the  positive  x-direction,  and  thus  k*r  = kx. 
In  Eq.  (1) 

k = m/c  ; c = (y/p)^'^^  , (2) 

s s 

and  u,  V and  w refer  to  the  displacements  in  the  x,  y and  z direc- 
tions respectively.  We  use  the  superscript  (O)  to  indicate  an  incident 

wave,  and  c its  wave  speed  in  the  matrix, 
s 

When  the  wave  impinges  on  the  scatterers,  part  of  the  incident  wave 
is  scattered  back  into  the  matrix  and  the  rest  is  refracted  into  the 
scatterer.  The  scattered  wave  will  also  be  of  SK-type.  The  displacement 

5 

component  of  the  scattered  wave  is  denoted  by  w = w while  that  correspond- 

F 

mg  to  the  refracted  wave  is  denoted  by  w = w . Since  the  incident  wave 


b . 
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has  time  dependence  given  "by  exp(-iu)t),  all  field  quantities  will  have 

the  same  time  factor  which  will  be  suppressed  henceforth.  The  dis- 
S F 

placement  field  w and  w satisfy  the  following  Helmholtz  equations 

(V^+k^)w®  = 0 ; (V^+k^)w^  =0  (3a,b) 

where  = 3^/3x^  + 3^/3y^,  and 

k^  = uj/c^  ; c^  = . (i*) 

The  index  f pertains  to  the  inclusion  material.  The  problem  at  hand 
reduces  to  computing  the  total  displacement  field  at  any  point  in  . 
the  matrix  satisfying  the  appropriate  boundary  conditions  on  the  sur- 
face of  the  inclusions  and  radiation  conditions  at  infinity. 

The  total  displacement  field  at  any  point  outside  the  scatterers 
is  the  sum  of  the  incident  displacement  field  and  the  fields  scattered 
by  all  the  scatterers.  This  may  be  written  as 

H 

w(r)  = w^(r)  + ^ w^(r-r^)  (5) 

i=l 

^ til 

where  is  the  field  scattered  by  the  i^'  scatterer  to  the 

point  cf  observation  r.  The  field  that  excites  or  impinges  on  the 

tfc  0 

1'  fcatterer  is  the  incident  displacement  w plus  the  field  scattered 

E 

* a-l  tfie  other  obstacles.  The  term  exciting  field,  w , is  used  to 

• -'-veer  the  field  actually  incident  on  a scatterer  and 

'.eikd  V produced  by  sources  at  infinity.  Thus, 
th 

-f  the  i ■ scatterer,  we  write 


i 

\ 

i 
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E/  N 
w^(r) 


0 x ' g 

w (r)  + y w.(r-r,)  ; a < |r-r. I < 2a 

j ~ ~J  _ 1 - -2. ' — 


vhere  'a'  is  the  radius  of  the  circle  circumscribing  a scatterer 


(see  Fig.  lb).  In  the  preceding  formulation,  ve  have  assimaed  that  the 


intercylindrical  spacing  is  such  as  to  avoid  any  overlapping  (no 


interpentration) . 


The  matrix  formulation  of  scattering  differs  from  the  eigen- 


function expansion  technique  in  that  the  same  basis  set,  for  our 


case  the  cylindral  wave  functions,  may  be  used  for  scatterers  of 


any  shape.  We  expand  the  exciting  field  and  the  scattered  field  in 


terms  of  basis  functions  iJ;  . auid  Fe  , 

n ' ’^n  ’ 


Wi(r)  = 


/ Re  il)  (r-r . ) 
Z_i  n n - -1 


j * n n ~ -j 


Where 


in<t> . 

(J)  = H [k Ir-r .Me  ^ 

n n ‘ - -1 ' 


in4>. 

’*'n  “ Jn[k|r-rj^|]e  ^ 


are  solutions  of  the  Helmholtz  equation,  and  a and  b are  undeter- 

n n 


mined  coefficients.  In  Eq.  (9)»  refers  to  the  angle  that 


makes  with  the  x-eixis  (see  Fig.  la)  and  Re  ili  refers  to  the  regular 

n 


part  of  i.e.,  the  part  that  is  regular  at  the  origin.  For 
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this  choice  of  the  "basis  set,  regular  part  of  is  different  from  the 

real  part  of 

It  has  been  shovra  that  if  the  total  field  outside  a scatterer  is  the 

s\un  of  the  incident  and  scattered  fields,  the  unknovn  scattered  field 

coefficients  can  be  related  to  the  incident  field  coefficients  throvigh 

17 

the  transition  matrix  (T-matrix)  . We  shall  extend  this  definition  to 

E S 

the  present  context.  Since  (w.  + v^)  is  the  total  field  in  the  matrix 

J J 


medium,  the  expansion  coefficients  of  the  field  scattered  by  the  j 


th 


scatterer  may  be  formally  related  to  the  coefficients  of  the  field  exciting 
the  J*”  scatterer  through  the  T-matrix: 


00 

b>3  = y tJ  aJ  . 

n nzn  m 

hf-od 


(10) 


The  elements  of  the  T-matrix  involve  surface  integrals,  which  can  be 
evaluated  in  closed  form  for  circtd.ar  and  spherical  geometry;  for 
obstacles  of  arbitrary  shape  they  can  only  be  evaluated  numerically. 

The  T-matrix  for  a single  scatterer  is  of  the  form 

T = -Re  Q(q"^)  (11) 

where  Re  Q and  Q are  matrices  which  are  functions  of  the  surface  S 

of  the  scatterer  and  of  the  natiire  of  the  boundary  conditions.  A detailed 

17 

derivation  of  Eqs.  (lO)  and  (11)  was  given  by  Waterman  • For  the  sake 
of  completeness,  expressions  for  Q-matrix  and  the  evaluation  of  the 
T-matrix  are  presented  in  the  Appendix  A. 

Substituting  Eqs.  (7)  and  (8)  in  Eq.  (6)  and  using  Eq.  (lO),  we  ob- 


tain 


C - 8 - 


Re  i|)  (r-r.)  = 

D > -i' 


I m 


It  can  be  seen  that  the  series  on  the  right  hand  side  of  Eq.  (12)  are 

til 

expressed  with  respect  to  the  center  of  the  j scatterer.  The  addi- 
tion theorem  will  be  employed  later  to  express  these  quantities  with 
respect  to  the  center  of  the  i^^  scatterer.  It  then  remains  to 
expand  the  incident  wave  also  in  the  form  of  a series  centered  at  the 
i"^^  scatterer.  By  employing  Fourier  series  expansion  in  complex  form. 


the  incident  wave  can  be  expanded  as 


xk  * X 

w^(r)  = y i^  Re  (r-r.)  e ^ 

~ n - '1 


With  Eq.  (13),  Eq.  (12)  can  be  written  as 


ik*r . y. 

Re  \Ij  (r-r.)  = e ^ y i^  Re  i|.'  (r-r.)  + 

n ~ ~i  p “ 'i 


j5^i  £ m 


T^  ip  (r-r.).  (II4) 
£m  £ - - j 


In  order  to  express  the  basis  function  iJ<j^(r-rj)  in  terms  of  V^^(r-r^), 
we  use  the  following  addition  theorem  (see  Fig.  2) 


,(k|r-r^|)e^  = X l^'^i  I ^ Vm^^iJ  ^ ® ^ ’^^5) 


After  some  simple  manipialations  and  using  Eq.  (9)»  ve  can  write  Eq.  (15)  as 


Substitution  of  Eq.  (I6)  in  Eq.  (ll()  yields 
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ik*r.« 

Re  (r-r,  ) = e ^ Re  \|;  (r-r . ) 


(17)  . 


N 


+ ^ (r.-r  ) Re  (r-r  ). 

Z-  Z_.Z_.Z^  to  m Jl-s  ~i  -j  S'  -1 

j#i  £ m s 


Multiplying  both  sides  of  Eq.  (1?)  by  Re  (r-r.)  and  using  the  orthogonality 

L **  1 

condition  on  the  angular  part  Eq.  (17)  can  be  reduced  to 


1 

a 

n 


ik*r . 

.n  ~ '1  ^ 
1 e + 


J’^i  £ m 


T"  a"^  (r  .-r  . ) , 
Jim  m £-n  'i  ~j 


(18) 


Thus,  by  employing  the  T-matrix,  ve  have  eliminated  the  unknown  scattering 

coefficients  b^,  resiilting  in  an  equation  for  a^  only.  In  view  of 

Eq.  (lO),  it  can  be  seen  that  the  scattered  field  coefficient  b^  can 

be  obtained  by  multiplying  Eq.  (l8)  by  the  transition  matrix  T.  The 

coefficients  a , or  b are  fvinctions  of  the  coordinates  r,,  6,  of  all 
n n J J 

the  scatterers,  although  this  is  not  written  explicitly. 

III.  STATISTICALLY  AVERAGED  WAVE  FIELDS 

We  now  wish  to  average  over  the  positions  of  all  the  scatterers. 

We  define  a probability  density  of  finding  the  first  scatterer  at  r^  and 
the  second  scatterer  at  r^  and  so  forth  by  p(r^,  r^,  Tj^).  The 

probability  density  may  be  written  as 

p(ri,  r^,  ....  Tjj)  = p(r^)  p(r^  r^,  r^jlr.) 


J-  5i> 


(19) 
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where  p(r^)  denotes  the  probability  density  of  finding  a scatterer 

at  r.  while  p(r.lr.)  denotes  the  conditional  probability  of  finding 

a scatterer  at  r,  if  a scatterer  is  known  to  be  at  r..  A prime  in 

the  first  of  Eq.  (19)  means  r^  is  absent  while  two  primes  in  the 

second  of  Eq.  (19)  means  both  r.  and  r are  absent. 

■"i  * j 

If  the  scatterers  are  randomly  distributed,  the  positions  of  all 
scatterers  are  equally  probable  within  volume  V,  and  hence  its  distri- 
bution is  uniform  with  probability  density 


p(r^) 


r . 

-1 


e V 


» 


r . 

'1 


i V . 


(20) 


In  addition,  the  assumption  that  the  cylinders  are  non  overlapping 
suggests  that 


P(r^|rj)  = 


1/V  , |r.-rj|  > 2a 


(21) 


where  'a'  is  the  radius  of  the  circumscribing  circle,  see  Fig.  1.  A 
suitable  correlation  in  the  position  of  cylinders  may  also  be  added  to 
(21).  The  correlation  is,  however,  omitted  here  for  simplicity. 

We  denote  the  conditional  expectation  of  a statistical  quantity 


f as 


f >i  = j"  ...  Jf  p(r^,  r^,  rjjjr.)  dr^  dr^  ...  dTj^ 


V V 


< f > 


/"•••  ^2’  *'*'»  •••'»  *^2’  ***'*  •••'» 


V V 


(22) 
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In  the  first  of  Eqs.  (22)  we  have  averaged  over  all  scatterers 
except  the  and  in  the  second,  over  all  except  the  i^^  and 

j , and  so  on. 

Multiplying  both  sides  of  Eq.  (l8)  by  the  probability  density 
given  by  (19)  and  using  Eqs.  (20),  (21)  and  (22),  we  obtain  the  con- 
figurational average  of  a^: 


1 

< a > . 
n 1 


ik*r . 
.n  - -1 
1 e 


tI  j < a.^  > ir  ) dr  (23) 

im  m £.-n  'J 


where  V denotes  the  volume  of  the  medi-um  excluding  the  volume  of 

K 

the  cylinder  of  radius  2a.  For  identical  scatterers,  Z can  be 
replaced  by  (K-1).  Equation  (21)  can  be  rewritten  as 


a .n 

< a >.  = 1 e 
n 1 


'i  ^ IW,  (-1)’^'^*'  T / < a*^  > ip  (r.,)  dr..  (2h) 

N ^ im  ' m ij  ^H-n  -j 


'^1 1 r ">  c 

' I »i  ■ ' ' / 

The  above  equation  indicates  that  the  conditional  average  with  one 


scatterer  fixed  is  given  in  terms  of  the  conditional  average  with  two 
scatterers  fixed.  This  presents  a difficulty  in  solving  for  the  un- 
known  coefficient  < a^  >.  However,  Lax  has  suggested  a quasicrystalline 
approximation , i . e . , 


< >4  > i J 

m xj  m j 


This  approximation  implies  that  there  is  no  correlation  between  the 
i^^  and  scatterers  other  than  that  there  shovild  be  no  inter- 
penetration. With  this  approximation,  Eq.  (,2h}  may  be  written  as 


i 
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1 .n 

< a >.  = 1 e 
n 1 


This  is  a system  of  integral  equations  for  the  unknown  coefficients 
< a^  >^.  Similar  expressions  for  < may  be  obtained  by  mxaltiply- 

ing  (26)  by  the  T-matrix. 


IV.  PROPAGATIOK  CHABACTERISTICS  OF  THE  AVERAGE  VAVES  IK  THE  MEDIUM 
We  now  try  a plane  wave  solution  for  the  average  < a > with  an 
effective  wave  number  K characterizing  the  composite  medivim: 


1 

< a >. 
m 1 


iK*r. 
.m  ' -1 

1 X e 
m 


The  averaged  wave  vector  K is  assumed  to  be  parallel  to  that  of  the 
incident  wave,  which  in  the  present  case  is  along  the  x-axis.  In  Eq. 
(27),  X^  is  an  unknown  constant.  Substituting  Eq.  (27)  into  Eq.  (26), 
we  obtain 


iK*r.  ik*r. 

.n  - -1  .n  - -1 

1 X e = 1 e + 

n 


¥22(-1)”*‘'^  i"x  / (r.Jdr. 

V ^ ^ £m  m J t-n  “ij  "j 


Equation  (28)  may  be  rewritten  as 


iKT,  ik*r 

.n„  * -i  .n 

i X e = a e 

n 


^*¥22<-i>'”*Vn).a/ 


The  integral  in  Eq.  (29)  is  a volume  integral  and  the  integrands 


iK-r. , 


and  i|)  (r..)  satisfy  the  wave  equation.  Thus  we  can  rewrite 
X -IJ 


the  integral  in  the  following  form 
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I. 


iK*r. 

e - ♦t(ry)dTj 


CD 

e?/ 


O iKT.,  p 

[V^{e  ' ^ ^ 


(30) 


I 


S -S^ 

00  2a 


3 


iKT..  2 

J 


vhere  Green's  theorem  has  been  applied.  refers  to  a circle  of  large 

radius  at  infinity  while  S refers  to  the  surface  of  the  circumscribed 
circle  of  the  scatterer. 

Using  the  following  plane  wave  expansion 


iK(r,cose  -r.cose. ) ^ i 

J J ^ ^ = y i"  J (Kr.J  e 

Z-.  n ij 


QS-OD 


the  surface  integral  on  can  be  evaluated: 


, - iKr.cose.  iKr  cose 


2a 


(31) 


iKr.cose  .t 

e ^ ^ o tr  [2k  J^(2Ka)H;(2ka)-2K  H.  (2ka)j;  (2Ka)  ] . 

k^-K^  t t t t 


(32) 


To  evaluate  the  surface  integral  on  S^,  we  replace 


asymptotic  value.  We  note  that  when 


-►  00 


- C ll*  - 


r .T. 
-1 


and  hence 


-*•  0 


J 


-ikr . T. 


ikr,  ite. 


iirkr, 


(33) 


(3M 


where  r is  the  tinit  vector.  Employing  Eqs.  (33)  and  (3^),  the  surface 
integral  on  S may  be  written  as 


r r . iKr . 


iKr.cose.  iKr,cose, 

J _ ^ J J _L 

3r, 


- e 


dS 


(35) 


2 2 
k -K 


J iTTkr. 

c J 


ar 


J 


ar,  ViTTkr 

t)  tJ 


In  Eq.  (29)  there  are  two  sets  of  terms  each  of  which  satisfies  a 
different  wave  equation  and  hence  both  can  be  equated  separately  to  zero. 
Thus  the  incident  field  term  and  the  surface  integral  at  infinity,  which 
satisfy  the  wave  equation  with  nvnnber  k cancel  each  other.  This  is  the 
so  called  extinction  theorem:  the  unperturbed  incident  wave  is  extinguished 
within  the  medium  by  waves  induced  at  the  boundary  of  the  system.  The 
remaining  terms  in  Eq.  (29)  with  Eq.  (32)  gives  the  relation  for  the  \in- 
known  conditional  averaged  amplitude: 


X 

n 


IZ 

£ m 


i““^X  T.  (JH.  ) 
m £m  £-n 


de . . 

J 


(36) 
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where  we  have  assumed  that  N and  V are  infinitely  large  with  H/V  = 
(finite),  and 

JH„  = 2kaJ.  (2Ka)H'.  (2ka)  - 2Ka  K.  (2ka)J'-  (2Ka)  . (37) 

t-n  Jt-n  t-n  t-n  *--n 


Similar  equation  can  also  be  obtained  from  the  average  scattered  field 
coefficient  ^ assuming  a plane  wave  solution  of  the  form 


iKr.cose. 

,1  .n  ^ 11 

< b >.  = 1 Y e 

n 1 n 


(38) 


where  Y is  an  unknown  constant, 
n 

In  Eq.  (36),  the  summation  integers  Z and  m varies  from  - <*>  to 
CO.  For  numerical  computations,  these  have  to  be  expressed  from  0 to 
This  can  be  achieved  by  formulating  certain  properties  of  the  T-matrix 
which  are  outlined  in  the  Appendix  B.  Using  these  properties  in  Eq.  (36), 
we  obtain 


eo 


TT-.  . 

/ T®  (X  +X  )(JH„^  +JH„  ) 

Z-  2 Jim  m -m  Jl+m  £-n' 


LJI=1 


.m 

+ T°  (X  -X  )(JH,  -JH.^„)  + — Htt  T®  (X  +X  „)JH 

Ji,m  m -m  Jl-n  Ji.+n  om  m -m  n 


CD 

£=1 

where  superscripts  ji,  and  o refer  to  even  and  odd  parts  of  T-matrix 
elements  (see  Appendix  B). 


(2) 


.Ji- 

~r  T®  X (JH„  +JK„^  ) + T X JH 
1/2  Ji^  o Jl-n  Jl+n  oo  o ni 


(39) 


Equation  (39)  is  a system  of  simultaneous  linear  homogeneous  equations 
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for  the  unknown  amplitudes  X^.  For  a nontrivial  solution,  we  require 

that  the  determinant  of  the  coefficients  vanishes  which  yields 

an  equation  for  the  effective  wave  number  K in  terms  of  k and 

the  T-matrix  of  the  scatterer.  This  is  the  dispersion  relation  for 

the  composite  medium.  It  can  be  shown  that  the  above  equation  is 

9 

similar  to  the  one  obtained  by  Bose  and  Mai  , for  circular  fibers, 
if  we  substitute  the  T-matrix  of  a circular  cylinder  given  in  the  Appen- 
dix A.  However,  Eq.  (39)  is  a general  expression  valid  for  any 
arbitrary  cross  section.  It  may  be  noted  that  2a,  the  distance  of 
closest  approach,  coincides  with  the  diameter  of  the  circ\alar  cylinders. 
The  T-matrix  is  the  only  factor  that  contains  infomation  about  the  exact 
shape  and  nature  of  the  cylinder.  If  the  elements  of  the  T-matrix  are 
known,  the  dispersion  relation  of  a composite  consisting  of  random 
distribution  of  obstacles  of  arbitrary  cross  section  can  be  obtained 
either  from  Eq.  (36)  or  Eq.  (39). 

V.  RAYLEIGH  LIMIT 

In  the  Rayleigh  limit  or  low  frequency  limit,  the  size  of  scatterers 
is  considered  to  be  small  when  compared  to  the  incident  wavelength.  It 
is  then  sufficient  to  take  only  the  lowest  order  coefficients  in  the  ex- 
pansion of  the  fields.  Neglecting  higher  order  terms  in  ka  and  Ka 
in  the  expansion  of  Bessel  and  Hankel  functions,  Eq.  (37)  may  be  reduced 
to 

JH  % |r[(f)"-2a^(k^-K^)£nka]  (39) 

li  7>o  il 

In  this  limit,  the  elements  of  T-matrix  are  known  in  closed  form  for 
various  simple  shapes  see  Ref.  [19],  and  they  are  given  here  for  the 
sake  of  completeness  for  circular  and  elliptical  cylinders: 
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LiMiL. 


Circular  cylinders 


T®o  = ^(d-l)k^a^[l-  ~ k^a^in  ka]  + 0(k^a’^) 


^00  = “ ■ *ln 


= 0 for  £,  n > 2 


Ellintical  cylinders 


TqO  ^ ^ J-(d-l)k^ab[l-  ^ k^ab  £n  ka]  + O(k^a^) 


I-L-rS.)  Lat^.  [ 1-  x ^ ^ ( a+b  )k^ab  £n  ka]  + 0(k^a^) 

8 a+mo  ^ i*  a+mb 


2^  (l-ra)(a+b)  ^^2ab[l-x  (a+b)k^ab  £n  ka]  + 0(k^a^) 

8 b+ma  _ k b+ma 


^00  = ° ’ ^£n 


= 0 for  £,  n > 2 


where  d = p^/p  and  m = v^/v  and  a and  b are  semi  major  and 
minor  axes  of  the  elliptic  cylinder,  respectively. 

In  view  of  Eqs.  (4o)  and  (i+1),  we  obtain  the  following  simultaneous 
linear  homogeneous  equations  from  (39)  for  the  tinknown  X^,  and 
X ^ for  thin  cylindrical  fibers 
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2-frn 
o 

2 2 
k -K^ 


2irn^ 

~2 


+ -jitT^(JH2+JH„)-T“^(JE 


J„-JH,))] 


’‘-l'  v2 
k 


2Trn  r ^ 

f l"n»” 


,*JH„)-M°^(JH^-JH„)) 


* -rlTn<JH2-JHo>-^ll<™2-^“ 


0*}  • 


(1*2) 


From  the  determinant  of  the  coefficients  X^,  and  X_^  ve  obtain 
the  following  dispersion  relations: 


Circular  cylinders 


where  c = Tra  n^ 


is  the  concentration  of  circular  cylinders  per  \init  area 


F 
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Elliptical  cylinders 


^-'^l  2UbTI)  1 -* 


2 2 

+ k a Jin  ka 


(-1 


^‘^l  (l-m)^(a+b)^ 


(mb+a) 
(l-ffi) (a+b] 
2(mb+a) 


-l+c^(d-l) 


+ O(k^a^) 


(a-D) 

(mb+a) 


(1*M 


where  c^^  = tab  n^  is  the  concentration  of  elliptical  cylinders  per 
unit  area. 

If  terms  of  order  k^a^JLn  ka  are  neglected,  it  is  seen  that  Eq.  (1*3) 

2h  7 

is  the  same  as  the  one  obtained  by  Hashin  and  Rosen  and  Bose  and  Mai  . 
For  the  same  order  of  approximation,  it  can  be  shown  that  Eq.  (1+1*)  agrees 

9 

with  the  results  obtained  by  Datta  for  low  concentration. 

The  value  of  K as  determined  by  Eqs.  (1*3)  and  (1*1*)  is  a real 

quantity  and  relates  to  phase  velocity  v^  of  the  composite  medium  given 

■'^p  = u)/K.  To  study  the  dependence  of  phase  velocity  on  concentration 

of  fibers,  we  consider  boron-aluminium  composite  with  density  ratio  d = 

2.53/2.72  and  shear  modulus  ratio  m = 25/3.87.  Using  these  values  in 

Eqs.  (1*3)  and  (1*1*),  we  computed  the  phase  velocity  v and  results  are 

P 

shown  in  Fig.  (3)  for  various  values  of  b/a  where  the  phase  velocity  v 

P 

is  normalized  with  respect  to  velocity  of  wave  propagation  c in  the 

s 

matrix.  The  general  tendency  of  the  phase  velocity  is  to  increase  with 
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concentration.  The  resiilts  also  indicate  that  as  b/a  decreases,  the 
phase  velocity  increases. 

2 2 

Defining  an  average  shear  modvilus  <ii>  = to  <.d>/K  where 
<p>  * cp^+(l-c)p  is  the  average  density,  ve  find  the  effective  shear 
modulus  of  the  composite  medium: 


<U> 

y 


l-c(l-m)/(l+m) 

l+c(l-m)/(l+m) 


(circular  cylinders) 


^ ^ l-c^(l-m) (a+b)/2{mb+a) 

-f-  " i+r^'(l-m)TaTbT7'2(mWa)  cylinders) 


(1*5) 


VI  DISPERSION  AT  HIGHER  FREQUENCIES 


To  study  the  response  of  the  composite  at  higher  frequencies, 
one  has  to  consider  higher  powers  of  ka  and  this  implies  that  a 
larger  number  of  terms  (X^)  must  be  kept  in  the  expansion  of  the  average 
field.  This  is  best  done  numerically.  Some  model  calculations  are 
presented  for  a al\iminium  matrix  reinforced  by  boron  fibers  of  elliptical 
cross  section.  The  material  properties  used  are  p = 2.72,  y = 3.87, 
p^  = 2.53,  and  = 25. 

For  values  of  ka  ranging  from  0.05  to  3.0,  the  determinant  of 
the  coefficients  of  X^  was  computed  nvunerically  retaining  13  simul- 
taneous, homogeneous  complex  equations  for  X g...XQ...X£.  The  elements 
of  the  T-matrix  were  computed  as  explained  in  Ref.  (17).  The  complex 
determinant  of  the  coefficients  was  calculated  using  the  Gaussian 
elimination  process  with  partial  pivoting.  For  a given  value  of  ka, 
the  roots  of  the  determinant  were  searched  in  the  complex  K-plane 
(K^+i  K^)  using  Muller's  method.  Good  initial  guesses  were  provided 
by  Eq.  (Itlt)  at  low  values  of  ka  and  these  could  be  used  systematically 
to  obtain  quick  convergence  of  roots  at  increasingly  higher  values  of 
ka.  The  real  part  determines  the  phase  velocity  v^  while  the 
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imaginary  part  determines  damping  or  attenuation.  The  results  of  these 
calculations  for  various  values  of  h/a  and  c are  ploted  in  Figs.  k-9. 

In  Figs,  k-6,  the  average  phase  velocity  in  the  composite  narmalized 
by  the  phase  velocity  in  the  matrix  material  is  plotted  versus  ka  for 
b/a  = O.Ji,  0.6,  0.8  and  1.0  and  for  c = 0.3,  0.5  and  0.7.  We  observe 
that  ■''■p/Cg  = k/K^.  For  b/a  «=  1.0,  these  results  agree  qixalitatively 
vith  those  obtained  by  Sutherland  and  Lingle^^  for  tungsten-aluminum 
composite.  Ting  and  Sachse  for  boron-epoxy  composites  and  Yew  and 
Jogi  for  steel-PD'J-i^  composites.  The  phase  velocity  increases  very 
slightly  at  low  frequencies  and  then  decreases  very  gradually  with 
increasing  frequency.  For  lower  values  of  b/a,  the  phase  velocity 
increases  much  more  rapidly  at  low  frequencies  and  then  tapers  off 
gradually  to  almost  the  same  values  as  for  b/a  =1.0.  It  may  be  noted 
that  our  results  indicate  that  the  phase  velocity  is  much  higher  for 
smaller  values  of  b/a  at  low  frequencies.  This  effect  is  more  pro- 
nounced at  higher  concentrations  where  the  velocity  also  drops  more 
sharply  with  increasing  frequency. 

In  Figs.  7-9*  we  have  plotted  the  coefficient  of  attenuation  a 
versus  ka.  We  chose  to  define  a as  o = itit  K^/K^  so  that  it  is 
dimensionless.  These  results  appear  to  be  new,  since  explicit  results 
for  the  attenuation  of  waves  in  composites  is  lacking  in  theoretical 
calculations  and  as  discussed  in  Refs.  15  and  16,  too  difficxilt  to 
measure  experimentally.  However,  the  results  presented  in  Figs.  7-9  agree 
qualitatively  with  those  predicted  by  Sutherland  and  Lingle^^.  The  at- 
tenuation increases  rapdily  at  first  with  increasing  frequency  (up  to 
ka  0.5)  and  then  decreases  to  almost  zero  at  high  frequencies.  There 
is  a tendency  for  the  attenuation  to  increase  rapidly  again  with  frequency 
This  can  be  observed  clearly  in  Fig.  7 for  c = 0.3.  In  these  results. 


p 


iMimuii.  ■ 


T 
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the  hehavior  is  qualitatively  the  same  for  different  values  of  h/a,  although 
the  increase  and  decrease  in  o becomes  more  pronounced  v’ith  decreasing  b/a. 

The  general  behavior  of  the  a versus  ka  curves  indicates  the 
possibility  of  existence  of  higher  modes  of  vibration  in  the  composite, 
or  what  is  referred  to  in  the  experimental  work  in  Ref.  [15]  as  the 
higher  pass  bands.  Although  Figs,  k-6  seem  to  indicate  that  there  is 
only  one  phase  velocity  in  the  composite  for  each  value  of  frequency, 
this  is  really  not  so  because  at  higher  frequencies,  there  exists  many 
values  of  K for  a given  ka.  These  multiple  roots  are,  however, 
difficiilt  to  obtain  in  our  root  searching  algorithm.  These  higher 
roots  simply  indicate  that  at  higher  frequencies,  the  effect  of  the 
individual  fibers  become  more  apparant  and  the  composite  behaves  less 
and  less  like  a continuum. 

The  frequency  at  which  a begins  to  decrease  is  usually  referred 

to  as  the  cut  off  frequency  of  the  first  pass  band.  It  must  be  noted 

that  the  attenuation  is  due  to  only  geometric  dispersion  or  scattering 

and  no  energy  is  dissipated  in  the  medium.  At  the  cut  off  frequency, 

a decreases  because  energy  begins  to  pass  into  the  second  pass  band. 

Our  results  seem  to  indicate  that  the  onset  of  the  second  pass  band  is 

q\aite  rapid  at  low  values  of  c as  indicated  by  the  rapid  decrease 

of  a to  0 for  b/a  = 1.0,  whereas  for  higher  values  of  c the 

15  22 

onset  is  more  prolonged.  Sutherland  and  Lingle  and  Moon  and  Mow 
offer  a similar  explanation  from  their  models.  It  should  be  mentioned 
here  that  o\ir  results  are  accurate  for  low  concentrations  because  of 
the  Lax's  quasicrystalline  approximation  employed  in  the  analysis 
which  neglects  the  correlation  between  the  fibers. 

There  are  a number  of  improvements  that  we  can  make  to  the  model 
composite  chosen  here.  The  reuidom  distribution  of  fibers  may  be 
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replaced  by  regular  arrays  of  fibers  since  this  is  true  in  practice,  and 

also  the  effect  of  correlation  between  fibers  must  be  taken  into  account 
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especially  at  high  concentrations.  As  observered  by  Peck  and  Sve  , 
the  effect  of  porosity  or  the  presence  of  voids  in  real  composites 
affects  the  attenuation  properties  quite  significantly  and  should  be  in- 


cluded in  oiu:  theoretical  model. 


APPENDIX.  A 


Transition  Matrix  for  a Single  Scatterer 


The  T-  or  transition  matrix  for  an  obstacle  bounded  by  the  surface 
S can  be  obtained  from  the  interior  and  exterior  Helmholtz  formulae 
for  the  scalar  vave  equation.  The  total  field  w outside  the  scatterer 


may  be  vritten  as 


v(r)  = v^(r)  + v^(r) 


(A-1) 


0 S 

vhere  v is  the  field  incident  on  the  scatterer  and  v is  the  scat- 
tered field.  We  then  have 


1“  j*  Ew’n'*7'g(r  (r’ )-g(r  |r’ )n’*P'v’' 


{-w  (r)  , r inside  S 

(A-2) 

w^(r)  , r outside  S 


where  g is  the  Green's  function  of  the  scalar  wave  equation  given  by 


(V^+k^)g  = -l+Tr6(r-r') 


(A-3) 


In  Eq.  (A-2),  w'  and  un'*V'w'  are  the  unknown  displacement  and  trac- 
tion respectively  on  the  siirface  S. 

We  expend  all  the  quantities  in  the  integral  formula  in  a set  of 
orthogonal  basis  functions  \i)^(r).  For  two  dimensional  scalar  waves, 
the  most  convenient  set  of  basis  functions  is  the  circular  wave  func- 
tions H^(kr)e^*^^.  Since  w^(r)  should  be  regixlar  for  r inside  S, 


we  have 


\r)  = y a Re  il>  (r)  = ^ a J (kr)e^ 

n n ~ n n 


(A-1*) 
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Also,  V should  be  outgoing  waves  at  infinity  and  be  regular  for  r 
outside  S.  We  then  have 

w^(r)  = 3^ b ijj  (r)  = X'b  H (kr) 

- Zj  n n - ^ n n 

n 


(A- 5) 


0 . 

In  the  above  expressions,  a^  are  known  constants  as  w is  given, 

but  b are  luiknowns. 
n 

The  expansion  of  the  Green's  functions  g in  terms  of  basis 
f\anctions  is 

r 


Re  i|)^(r')  , |r|  > jr' 


g(rlr')  = 


n 


lir 


^ii/^(r')  Re  , |rj  < |r' | . 


(A-6) 


Substituting  Eqs.(A-l+),  (A-5)  and  (A-6)  into  (A-2)  and  using  the  ortho- 
gonality of  the  circular  wave  functions , we  obtain 


J [v'n'  -V'  ' )-il'^(r ' ) n'*V'w']dS 


(A-7) 


and 


b = 7"  r [w'n'*V'  Re  (r' )-Re  iji  (r')  n'*V'w']dS.  (A-6) 

n 4 y n - n ~ 


For  an  elastic  inclusion  with  shear  modulus  the  relevant 

boundary  conditions  are 


w _ = w _ „ 

•S  f 'S 


(A-9) 


and 


vn*Vw|  = ;j.n*Vw 


f-  f'S 


(A-IO) 
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The  unknown  field  inside  the  inclusion,  can  be  expanded  as 


V =3^0  J (k»r)  (A-ll) 

f n n f 

n 

where  k^  pertains  to  the  wavenumber  for  the  inclusion  material.  The 
traction  is  given  by 


n'  ‘-Vw^ 


n 


, r inside  S.  (A-12) 


Substituting  Eqs.  (A-ll)  and  (A-12)  in  Eqs . (A-7)  and  (A-6)  using  the 
boundary  conditions  (A-9)  and  (A-10),  we  can  write 


a 

n 


(A-13) 


m 


(A-ll*) 


where  we  have  defined 

Snn=  I r(({i)n(r)  • [Reiii^(k^r  )7i|/^(kr  )-(p^/p ) VReil)^(k^r  )\|/^(kr ) ] ] V[r-r(((>)  ] jdi}) 
® (A-15) 

Equations  (A-13)  and  (A-llt)  can  be  solved  to  yield,  in  matrix  notation, 

b = -(Re  §)9*^a  = -T  a.  (A-l6) 


Thus,  the  T-matrix  of  transition  matrix  for  an  obstacle  relates  the  coef 
ficients  of  the  field  incident  on  a scatterer  to  the  coefficient  of 
the  field  scattered  by  it  and  is  given  by 

T * -Re  Q (A-17) 


and 


n(r) 


For.  a circular  cylinderical  scatterer  of  radius  a,  r(4>) 
= r.  Equation  (A-15)  then  gives 


a 


= -(wa/l»)[J^(k^a)k  H;(ka)-(v^/p)k^ 


(A-18) 


Substitution  of  Eq.  (A-18)  into  Eq.  (A-IT)  yields  the  T-matrix  for  a circular 
elastic  cylinder 


vJ^(k.pa)kJ^Jka)-vi^J^(ka)k^J^(k^a) 
^mn  ~ vi^^(ka)k^J^(k^a)-uJjj(k^a)kH^(ka) 


{A-19) 


For  a scatterer  of  arbitrary  shape,  can  be  determined  from  (A-15) 


by  numer' 


integration. 


APPSrPIX  E 


Relations  r-'-,--.'-,'  e-»irer.tt  aeflned  in  tvo  different  basis  sets 

In  fornralatiiif  tx.*_  r scattering  problen:,  we  have  chosen  the 

complex  ang‘wilar  funct-cns  construct  the  basis  functions 

in  favor  of  the  real  a:.gu_ar  functions  cos  n(f  smo  sin  nip.  The  reason 
for  this  choice  is  that  the  addition  theorem  for  cylinder  functions  in 
the  form  of  Eq.  (15)  is  more  convenient  and  less  awkward  than  that  in 
terms  of  the  sines  and  cosines.  To  make  such  a choice,  we  have  paid 
the  price  that  the  real  part  of  the  basis  set  cannot  be  constructed 

by  simply  replacing  by  J^.  Had  the  real  angular  functions  been 

used,  "real  ip^"  and  "regular  ipj^"  are  then  equal. 

In  numerical  computation  of  the  T-matrix  elements,  it  is  convenient 
to  make  "real"  synonjrmous  with  "regular"  for  reasons  that  we  cannot  go 
into  here.  Thus  it  becomes  necessary  to  compute  the  T-matrix  in  a 
basis  set  with  real  angular  fxinctions  and  then  convert  it  into  a 
basis  set  with  complex  angular  fvinctions. 

To  distinguish  between  the  Q-  and  T-matrices  defined  in  the  two 
basis  sets,  we  shall  use  additional  superscripts  o,v  = e,  0 to 

denote  the  even  and  odd  angular  functions. 

0 S 

We  expand  w and  w in  the  two  basis  sets  and  obtain 


00 
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where 


a°  = itan-(-l)  a_j]//Cj 


(B-2) 


In  Eq.  (B-1),  is  the  Neumann  factor  = 1 and  = 2 for  £ > 0; 
i - 0,  1,  2,  ...  for  a = e and  Z = 1,  2,  3,  ...  for  o = o. 
Similarly,  we  write 


w^(r)  = ^ b H (kr)  e^' 

- z_)  n n 


‘I 


cos  £<ti  ; o = e 


sin  £.({i  ; o = o 


(B-3) 


The  coefficients  b^^  and  b^^^  are  related  as  given  by  Eq.  (B-2). 
Proceeding  as  in  Ref.  17,  ve  obtain 


3 = I T a 

n nm  m 

in=— “ 


(B-M 


b'^  = Z Z a'* 

n nm  m 

m V 


where 


(B-5) 


(B-6) 


C = -II 


(B.T) 


Let  the  surface  of  a scatterer  be  expressed  by 
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r-r(4i)  = 0 on  S 

The  corresponding  Q-  matrix  elements  are 


(B-8) 


2tt  ^ y 

Q f /[J  (k^r)  K'(kr)-  — H (kr) 

nm  k J I m f n pn  ml 


J (k^r)  H (k; 

^2(^)  a<p  v m . n 


.r)jei(”'-^)'^i 


{B-9) 


2^  . p 

qOv  = f /e  /e  r [j  (k^r)  H'(kr)-  — K„(kr)  J'{k^r)lF 
^nm  h n mji  m f'n  pn  n:-‘  nm 


1 d 


r^(*)  "♦ 


— r($)  H^(kr)[n  R_-  m ) r((^)d<;> 


where 


and 


nm 


nm 


nm 


nm  p nm 


'} 


cos 

n(}i  cos  m(fi 

cos 

nip 

sin  m(J>  \ 

sin 

ncfi  cos.miji 

sin 

nip 

sin  mip  / 

-sin  n<)> 

cos  m4i 

sin 

nip 

sin  m4>  \ 

cos 

n<t> 

cos  m<f 

cos 

nip 

sin  m4'  / 

-cos  n(ti 

sin  m(}) 

cos  n4i  cos  m4i  \ 

‘-sin  n<f‘ 

sin  m(t) 

sin  nif)  cos  m(})  / 

(B-10) 


(B-11) 


(B-12) 


(B-13) 


In  Ec?. (B-9)  and  (B-IO),  the  prime  denotes  derivative  with  respect  to 
r,  and  k,  and  k are  the  wave  numbers  in  the  inclusion  and  matrix 


materials  respectively. 

From  Eqs.  (B-9)  and  (B-IO)  for  boundaries  with  the  following 
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symmetry 


r(e)  = r(2Tr-e)  = r(Tr±e)  , 


(B-lM 


we  obtain 


nm 


Q = 0;o7^v;n  + m odd 

pTn 


(B-51) 


Defining 


.oo 


Q H Q 
nm  nm 


we  have 


nm 


2 nm  -n,m 


(b-i6) 


Q, 


nm 


7^{Q  -(-ifQ  ^ 

Yt  nm  -njm 


(B-17) 


(Q  ) 


/T  /r 


nm 


2 run  -n  >in 


■{B-i8) 


ana 


run  rt  nm  -n,m 

n 


(B-19) 


The  symmetry  implied  by  Eq,.  (B-lU)  is  applicable  to  all  quadratic  surfaces, 
for  example,  the  fibers  of  elliptic  cross  section  considered  in  the 
numerical  calculations.  Using  Eqs.  (B-16)  - (B-19)»  the  relation  between 
the  T-matrices  of  the  two  basis  sets  is  given  by 


T 

n,-m 


(B-20) 


(E-21) 


In  arriving  at  Eqs.  (B-20)  and  (B-21)  the  symmetry  property  of 
the  T-matrix  has  been  invoked.  Equations  {B-20)  and  (B-21)  are  used  to 
change  the  s\mmiations  in  Eq.  (36)  from  -“>  to  +“  to  positive  values 
only  as  given  by  Eq.  (39)- 

It  may  be  pointed  out  that  the  analytical  properties  presented  in 
Eq.  (B--.6)  - (B-21)  are  not  confined  to  the  scalar  wave  equation  but 
are  equally  applicable  to  the  vector  wave  equation  and  the  elastic 
wave  equation  since  the  definition  of  the  T-matrix  remains  the  same. 

ACKMO\-?LEDGEI.IEIfrS 

This  research  was  supported  by  the  National  Science  Foundation 
through  the  Materials  Science  Center,  Cornell  University,  and  by  the 
Air  Force  Office  of  Scientific  Research. 


1.  Lord  Rayleigh,  Phil.  Mag.  1*7 » 375  (1899) 

2.  L.L.  Foldy,  Phys.  Revs.  67.  107  (l9l*5) 

3.  M.  Lax,  Rev.  Mod.  Phys.  287  (1951) 

1*.  M.  Lax,  Phys.  Revs.  621  (1952) 

5.  V.  Twersky,  J.  /‘coust.  Soc.  Am.  29.  209  (1957) 

6.  V.  Twersky,  J.  Hath.  Phys.  3,  700  (1962) 

7.  P.C.  VJaternan  and  R.  Truell,  J.  Math.  Phys.  £,  512  (I96I) 

8.  P.C.  Karl,  Jr.  and  J.B.  Keller,  J.  Math.  Phys.  537  (I96I) 

9.  S.K.  Bose  and  A.K.  Mai,  Int.  J.  Solids  and  Struct.  £,  1075  (1973) 

10.  A.K.  Mai  and  A.K.  Chatterjee,  J.  Appl.  Mech.  11,  6l  (1977) 

11.  S.K.  Datta,  J.  Apnl.  Mech.  L2,  I65  (1975) 

12.  J.D.  Achenbach,  "Waves  and  Vibrations  in  Directionally  Reinforced 
Composites",  in  Comnosite  Materi£.ls  (L.J.  Broutman  and  R.H. 

Krock,  eds.)  Volume  2,  Academic  Press,  New  York  (1972) 

13.  J.C.  Peck,  "Stress  vave  propagation  and  fracture  on  composites", 
Symp.  Dynam.  Compos.  Mater.  ASKE  (1972) 

ll*.  F.C.  Moon,  "Wave  propagation  and  impact  in  composite  materials", 
in  Composite  Materials  (L.J.  Broutman  and  R.R.  Krock,  eds.) 

Volxme  75  Academic  Press,  Kew  York  (1975) 

15.  H.J.  Sutherland  and  R.  Lingle,  J.  Comp.  Mater.  L90  (1972) 

16.  T.R.  Tauchert  and  A.K.  Guzelsu,  J.  Appl.  Mech.  39,  98  (1972) 

17.  P.C.  Waterman,  J.  Acoust.  Soc.  Am.  i<5,  lUl7  (1969) 

18.  B.  Peterson  and  S.  Strom,  J.  Acoust.  Soc.  Am.  56,  771  (197^) 

19.  V.V.  Varadan  and  V.K.  Varadan,  Low  frequency  expansions  for  the 

scattering  of  scalar  waves  by  cylinders  of  simple  shapes.  Materials 
Science  Center  Report  ^ 2937  , Cornell  University,  Ithaca, 

Hew  York  (1977) 

20.  C.S.  Ting  and  W.H.  Sachse,  private  communications. 


C.H.  Yew  and  P.K.  Jogi,  Proceedings  of  the  12th  Annual  meeting  of 
Society  of  Engineering  Science,  University  of  Texas,  Austin,  Texas  ,1*97 
(1975) 

- C - 33  - 


FIGURE  CAPTIOKS 


Fig. 

1. 

Fig. 

2,. 

Fig. 

3. 

Fig. 

1*. 

Fig. 

5. 

Fig. 

6. 

Fig. 

7. 

Fig. 

8. 

Fig. 

9. 

Geometry  of  parallel  cylinders  and  circumscriting  circle 
Translation  of  coordinate  system  of  the  and  i^^  cylinders 
Normalized  phase  velocity  vs  concentration  - Rayleigh  limit 
Normalized  phase  velocity  vs  ka  for  c = 0.3 
Normalized  phase  velocity  vs  ka  for  c = 0.5 
Normalized  phase  velocity  vs  ka  for  c = 0.7 
Attenuation  coefficient  vs  ka  for  c = 0.3 
Attenuation  coefficient  vs  ka  for  c = 0.5 
Attenuation  coefficient  vs  ka  for  c = 0.7 


b/a  = O.U 


Fig.  ^ ftlRMALIZED 


o 


Ahenuation 


Attenuation 


Attenuation 


